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ABSTRACT 



Physical dust models are presented for 65 galaxies in the SINGS survey that are 
Q^ | strongly detected in the four IRAC bands and three MIPS bands. For each galaxy 

we estimate (1) the total dust mass, (2) the fraction of the dust mass contributed by 
■ PAHs, and (3) the intensity of the starlight heating the dust grains. We find that 

spiral galaxies have dust properties resembling the dust in the local region of the Milky 
£N| . Way, with similar dust-to-gas ratio, and similar PAH abundance. The observed SEDs, 

including galaxies with SCUBA photometry, can be reproduced by dust models that do 
! not require "cold" (T < 10 K) dust. 

The dust-to-gas ratio is observed to be dependent on metallicity. In the interstel- 
lar media of galaxies with Aq = 12 + log 10 (O/H) > 8.1, grains contain a substantial 
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fraction of interstellar Mg, Si and Fe. Galaxies with Ao < 8.1 and extended H I en- 
velopes in some cases appear to have global dust-to-gas ratios that are low for their 
measured oxygen abundance, but the dust-to-gas ratio in the regions where infrared 
emission is detected generally appears to be consistent with a substantial fraction of 
interstellar Mg, Si, and Fe being contained in dust. The PAH index qpah ~~ the frac- 
tion of the dust mass in the form of PAHs - correlates with metallicity. The nine 
galaxies in our sample with Ao < 8.1 have a median qpah = 1-0%, whereas galaxies 
with Ao > 8.1 have a median qpAn = 3.55%. The derived dust masses favor a value 
A^co ~ 4 x 10 20 cm _2 ( K km s" 1 ) -1 for the CO to H2 conversion factor. Except for some 
starbursting systems (Mrk 33, Tolo 89, NGC 3049), dust in the diffuse ISM dominates 
the IR power. 

Subject headings: ISM: dust, extinction — ISM: general — galaxies: abundances - 
galaxies: general — galaxies: ISM — infrared: galaxies 



Introduction 



Interstellar dust affects the appearance of galaxies, by attenuating short wavelength radiation 
from stars and ionized gas, and contributi ng infrared (IR ), far-infrared (FIR), submm, mm, and 
microwave emission [for a recent review, see iDraind (|2003l )]. Dust also is an important agent in the 
fluid dynamics, chemistry, heating, cooling, and even ionization balance in some interstellar regions, 
with a major role in the process of star formation. Despite the importance of dust, determination 
of the physical properties of interstellar dust grains has been a challenging task - even the overall 
amount of dust in other galaxies has often been very uncertain. 



The SINGS galaxy sample (jKennicutt et alj 12003 ) comprises 75 galaxies with a wide range 
of morphological types, all at distances allowing the galaxies to be resolv ed even at 160//m. Ob- 
servations of this sample with all instruments on Spitzer Space Telescope (jWerner et al.l 12004! ) are 
providing a new perspective on the dust abundances and dust properties in a diverse set of galaxies. 

For 67 of t he SINGS ga l axies we have positive detections in the 4 bands of the Infrared Array 
Camera (IRAC .lFazio et a l. 2004)) and the three bands of the Multiband Imaging Photometer for 



Spitzer (MIPS; iRieke et al.l (|2004! )). spanning wavelengths from 3.6/mi to 160^m, the wavelength 



range within which dust radiates most of its power, including the PAH emission features that are 
prominent in many galaxies; 65 of these galaxies are used in the present study. For 11 galaxies we 
have global fluxes measured with the IRS 16^m peak-up filter. For 17 of the 65 galaxies, global 
850/iin flux measurements are also available from the SCUBA camera on the JCMT . 

The IR emission from dust depends not only on the amount of dust present, but also on the 
rate at which it is heated by starlight. We develop a technique for estimation of dust masses and 
starlight intensities in galaxies, and apply it to estimate the total dust mass M<-j ust in each of the 65 
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galaxies in this study. For 61 galaxies we are also able to estimate a "PAH index" qpAn ~~ denned 
here to be the mass fraction of the dust that is contributed by polycyclic aromatic hydrocarbon 
(PAH) particles containing < 10 3 C atoms. In addition, we obtain information on the starlight 
intensity distribution - the mean starlight intensity scale factor (U), and the fraction of the dust 
heating that takes place in regions of high starlight intensities, including photodissociation regions. 

The paper is organized as follows. The sample selection and observational data are summarized 
in <J21 the physical dust models that are used are described in <|3j and the model-fitting procedure 
is given in In order to validate the model-fitting procedure, we first limit attention to the 
subset of SINGS galaxies for which SCUBA data are available (we will refer to these as the SINGS- 
SCUBA galaxies). In ^5] we demonstrate that the dust modeling procedure works well for most 
of the SINGS-SCUBA galaxies, and we present dust masses, dust-to-gas mass ratios, PAH index 
values, and starlight intensity parameters obtained from these model fits. For the SINGS-SCUBA 
galaxies, ou r models do not requ ire any contribution from the very cold (< 10 K) dust that has been 



reported bv iKrugel et a 
dust that 



Galliano et al. 



that Dumke et al 



(|2003l . 



(1199^ for the spiral galaxies NGC 6156 and NGC 6918, or the 5 - 9 K 



20051 ) reported for low metallicity dwarf galaxies, or the 4-6 K dust 



(J2004J) reported for the edge-on-spiral NGC 4631. The global SEDs do not show 



eviden ce for the submillimeter excess seen in the faint outer regions of NGC 4631 (jBendo et al 



2006bJ). 



Submillimeter measurements are valuable to help constrain the dust models, but are unfor- 
tunately not available for most galaxies. In $6] we use the SINGS-SCUBA galaxies to develop a 
"restricted" fitting procedure that can be applied to galaxies for which submillimeter fluxes are un- 
available. For the 17 galaxies with global 850/im fluxes from SCUBA, we show that the restricted 
fitting procedure, without employing SCUBA data, recovers dust masses and starlight intensities 
that are usually within a factor of two of the values obtained when the SCUBA data are used. 

The restricted fitting technique is then applied to the remaining 48 galaxies to estimate dust 
masses and starlight intensities. Dust-to-gas mass ratios for the overall sample are presented in 
§7-H where it is found that some of the SINGS galaxies have very low dust-to-gas ratios, and 
that the dust-to-gas ratios are correlated with both gas-phase metallicity and galaxy type. The 
dust-to-gas ratios depend on the assumed value of Xqo, the ratio of H2 column density to CO 
line intensity; in £19.41 we argue that the distribution of derived dust-to-gas ratios supports Xqo ~ 
4 x 10 20 cm- 2 (Kkm a -1 ) -1 . 

Values of the PAH index (/pah are estimated for 61 galaxies in §7.41 Many galaxies are found 
to have low values of qpAK, and qpah is found to be strongly correlated with gas-phase metallicity. 
The starlight parameters for the 65 galaxy sample are presented and discussed in §[BJ The results 
are discussed in $9j and summarized in §101 
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2. Data 



2.1. IRAC and MIPS 



Spitzer Space Telescope w as used to observe all 75 galaxies in the SINGS sample, following the 
observing strategy described in lKennicutt et al,l (|2003l ). The observations in the 4 IRAC bands (3.6, 
4.5, 5.7, 7.9f im) and 3 MIPS bands (24, 71, 160/mi) have been reduced following th e procedures 



descri bed in iDale et al.l (|2005l ). The extended source aperture corrections found in iReach et al 



(|2005l ) were applied to the IR AC data. Global IRAC and MIPS photometry for these galaxies is 
taken from [pale et al.l ( 2007 ). which included (usually modest) corrections for aperture size and 
nonlinearity of the 71 /im detectors. 



Ten of the SINGS galaxies were excluded from the present study: 



NGC 584 (E4) may be contaminated by a background source, and also has low surface 
brightness in the MIPS bands compared to the foreground cirrus. 

NGC 3034 = M 82 (10, starburst): Saturation of the bright core of this galaxy precludes 
global photometry. 

NGC 1404 (El): An off-center IR source contributes a significant fraction of the MIPS flux, 
but it is not known whether or not it is a foreground/background object. 

NGC 4552 (E0) has low surface brightness compared to the foreground cirrus emission, and 
the MIPS 160/im photometry is very uncertain. 

The dwarf galaxies M81dwA and Ho IX are not detected at A > 5.8 /im. 

The dwarf galaxies DDO 154 and DDO 165 are only marginally (< 3<r) detected at 71 and 
160/iin, and M81dwB is a < 2a detection at 160/iin. 

The dwarf starburst galaxy NGC 1377 shows strong silicate absorption at 9.8/xm and 18/tm 
(jRoussel et al.ll2006l ). The objective of the present study is to model the IR spectra using a 
simple model that assumes the galaxy to be optically-thin at A > 3.5 /jm; because NGC 1377 
is clearly not optically thin, it has been removed from the sample. 



The sample studied here therefore consists of 65 galaxies. MIPS 71/im photometry for NGC 7552 
is not usable because of detector saturation, but NGC 7552 is retained in the sample because it is 
well-observed at other wavelengths and with other instruments. 

Most of the dust mass in galaxies is relatively cool, radiating primarily at wavelengths in the 
MIPS 71 fim and 160 /im bands. Dust grains that are warmer, either because they are in regions 
with intense starlight, or because they are very small grains immediately following absorption of a 
stellar photon, may radiate strongly in the MIPS 24 /im band and shortward. The variety of dust 
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Fig. 1. — MIPS color-color plot for the 65 galaxies. The blue shaded area shows the region sampled by the 17 
galaxies for which global SCUBA fluxes are known (plotted as blue solid symbols) - these only partially cover the 
color-color space populated by the SINGS sample (green shaded area). Galaxies without SCUBA fluxes are shown 
in red. Galaxies with relatively extreme flux ratios are labelled. 



conditions found in the SINGS sample is evident from the distribution of the 64 galaxies in MIPS 
color-color space, shown in Figured) {vF^ji^J ^(^i ? i/)i60 j um varies by over a factor of 10 within 
our sample, and (vF u )24 fl m/{i / F v )ji fim varies by a factor of 8. These variations in MIPS colors 
presumably reflect galaxy-to-galaxy differences in the intensity of the starlight heating the dust 
grains, although galaxy-to-galaxy differences in dust composition may also play a role. It is evident 
from the scatter in Figure[T]that the MIPS 24/71 color is not simply a function of the 71/160 color: 
the MIPS colors cannot be described by a one-parameter sequence of, e.g., "dust temperature". 
Realistic dust models will therefore need to vary at least two parameters to account for the range 
in observed MIPS colors. 

The SINGS sample does not include any powerful AGN, and in all SINGS galaxies the global 
flux in the IRAC 3.6^m band is starlight-dominated. Neglecting extinction, we can estimate the 
"nonstellar" flux in longer-wavelength bands by subtracting the estimated stellar contribution, 
assuming a 5000 K blackbody. For IRAC band 4, the IRS 16/im Peakup band (see ^2.21 below) and 
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the MIPS 24/mi band, the nonstellar fluxes are: 
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Th e coefficients 0.260 and 0.0326 in eq. (pQ) and ([3]) are close to the values 0.232 and 0.032 estimated 
bv lHelou et al.l (|2004l ). Some stars (e.g., cool red giants, and AGB stars) have circumstellar dust 
that is in many cases warm enough to radiate strongly at wavelengths from 8-30/im. If such stars 
are sufficiently numerous, the result will be a "stellar" contribution to the 8, 16, and 24/mi fluxes 
that will be larger than assumed in eq. (JTHSl). In the absence of extinction at 3.6/mi, equations 
(Q~H3|) will therefore overestimate the flux contributed by interstellar dust. 



2.2. IRS 16^111 

In addition to IRAC and MIPS photometry, nine of the SINGS galaxies have been imaged 
using the IRS "blue peak-up" detector array, with nominal wavelength A = 16 /.mi. The IRS blue 
peak-up observations were taken as part of the IRS instrument calibration program and consist of 
a grid of peakup images covering the full extent of each galaxy. The reduced images, processed 
by IRS pipeline version S14, were mosaicked using the MOPEX program (ver. 030106) to remove 
cosmic rays after subtraction of a custom dark image. This custom dark image was created using 
images taken beyond the visible extent of the galaxy. In some cases, the best dark subtraction was 
found to require allowing the dark image to vary with time. "Sky" (zodiacal light and Galactic 
cirrus) contributions were also removed by this custom dark subtraction ste p. The total 16 u m. flux 



was measured using the same methods used for the IRAC and MIPS fluxes ()Dale et al.ll2005l ). The 
measured global fluxes were then scaled down by ~15%, to account for a recent SSC calibration 
update (~10%), and the appropriate aperture correction factor (~5%). Fluxes for these 11 galaxies 
at and 24/mi are given in Table [TJ The flux density ratio {F u )iq/ {F u )24 varies considerably 

among these galaxies, from a low of 0.376 ± 0.033 for NGC 7552 to a high of 1.21 ± 0.12 for 
NGC 3190. 

Figured shows P w = {vF^) 16 /((uF v ) 71 + {uF v ) leo ) vs. P 24 = {vF™) 2 ±/{(vF v ) n + (vF v ) 160 ) 
for the ten galaxies for which we have 16, 24, 71, and 160/um photometry, where F™ s is the nonstellar 
contribution from eq. ([HE]). The flux ratio P\q remains essentially constant at P\q ~ 0.112 ± 0.017 
while P24 varies by almost a factor of 3, from 0.06 to 0.18. This suggests that the 16p,m flux 



is not emitted by the same grains as the 24/mi flux. In the dust model of iDraine k, Lil (|2007l . 
hereafter DL07) , a minimum level of 24/im emission is produced by single-photon heating of PAHs 
or ultrasmall grains (with (^L^^/Ltir ~ 0.03 — 0.06, depending on gpAH _ see Fig. 15 of DL07), 
but P24 increases as the starlight intensity increases. Figure [2^, suggests that the 16/xm emission 
may be dominated by single-photon heating of PAHs, remaining relatively constant even as P24 
rises due to increased starlight heating. This is supported by Figure [2)3, which shows that the ratio 
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Fig. 2.— (a) Pi 6 = (isF™) lfi /((vF v } 71 + (uF v ) 160 ) versus P 24 = (isF™} 24 /((vF v } T1 + (uF v ) 1(i0 ) for the 10 galaxies 
imaged with the IRS 16/xm peak-up filter for which we also have global MIPS fluxes (NGC 7552 is excluded), (b) 
{vF^ s )7.g/ {vF„ s )i6 vs P24, Both the ratio of nonstellar 16/xm emission to total IR power and the ratio of nonstellar 
7.9/xm emission to nonstellar 16/xm emission remain constant to within ~ 15%, even as P24 varies by a factor of 3. 



of 7.9/xm emission to 16/xm emission is nearly constant for this sample, varying by only ~ ±15% 
over the eight galaxies, consistent with the idea that bo th are primarily the resu l t of s ingle-photon 



heating of PAHs. This is consistent with the finding bv lForster Schreiber et al.l (|2004l ) that, when 
averaged over disks of spiral galaxies, the 15^m continuum intensity is proportional to the PAH 
intensity in the 5-8/Lim region. 

Table 1. Global Fluxes (Jy) in IRS 16/im Peak-Up Band 



galaxy 




(F v ) 24 






Jy 


Jy 




NGC 0337 


0.375 ± 0.029 


0.678 ± 0.028 


0.553 ±0.049 


NGC 3190 


0.322 ± 0.030 


0.267 ±0.011 


1.21 ±0.12 


NGC 3351 


1.39 ± 0.23 


2.583 ±0.12 


0.538 ±0.091 


NGC 3521 


4.86 ± 0.59 


5.51 ±0.22 


0.88 ±0.11 


NGC 3627 


4.92 ± 0.58 


7.42 ±0.30 


0.663 ± 0.082 


NGC 3938 


0.849 ± 0.068 


1.087 ±0.045 


0.781 ±0.070 


NGC 4321 


2.64 ±0.43 


3.34 ±0.13 


0.79 ±0.13 


NGC 4536 


1.74 ±0.14 


3.46 ±0.14 


0.503 ±0.044 


NGC 5055 


5.09 ± 0.43 


5.73 ±0.23 


0.888 ± 0.083 


NGC 6946 


12.8 ± 1.02 


20.37 ±0.81 


0.628 ± 0.056 


NGC 7552 


4.01 ± 0.31 


10.66 ±0.44 


0.376 ± .033 
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2.3. SCUBA 

The dust properties in a galaxy will be best determined if submillimeter photometry is available 
to constrain the quantity of cool dust present. SCUBA observations at 450//m and/or 850/mi are 
available for 26 of the 65 galaxies but not all of these are usable in the present study: 

1. NGC 4254, NGC 4579, NGC 5194, and NGC 6946 were scan-mapped, and the data processing 
removes an unknown amount of extended emission. Thus global 850/im (or 450/mi) fluxes 
are not available for these galaxies. We restrict ourselves to data obtained in "jiggle map" 
mode, except for the scan-map observations of NGC 5195, for which background subtraction 
is thought to be more reliable because of the smaller angular extent of NGC 5195. 

2. For some galaxies the SCUB A map did not fully cover the region where 160/mi emission is 



detected. iDale et al.l (|2005l ) multiplied the observed SCUBA fluxes by a correction factor 
based on the fraction of the 160/mi emission coming from the region mapped by SCUBA. 
Because the 850/xm/160/xm flux ratio may be spatially variable, we choose here to not use 
the 850/im flux estimates for galaxies where the corrrection factor exceeds 1.6 (NGC 1097, 
NGC 4321, NGC 4736, and NGC 5033). 



3. For th e Sombrero galaxy, NGC 4594, the 850/mi flux is dominated by the AGN (jBendo et al 



2006al ). and the extranuclear 850/im flux is not reliably measured, so we cannot use the 



SCUBA data to constrain the dust model for this galaxy. 

We thus have 17 galaxies where the dust models are constrained by 850/im SCUBA data (with 
450/mi fluxes for 3 of the 17) in addition to IR AC and MIPS da ta. The global fluxes -F bs,ft an d 



1-a uncertainties for each band b are taken from lDale et al.l (|2007l ). 



The 160 galaxies for which we have global SCUBA and MIPS photometry are indicated in 
Figure [H While the galaxies with SCUBA photometry span a considerable portion of the populated 
region in MIPS color-color space, a number of the SINGS galaxies have MIPS colors that are outside 
of the range seen among the galaxies for which SCUBA data are available. For example, NGC 4725 
or NGC 5408, with (vF v ) 71^/(^^)160^111 = 0.33 and 3.1, respectively, fall outside the range 0.6- 
2.4 observed for this ratio among the SINGS-SCUBA galaxies. The dust models discussed below 
are able to reproduce the IRAC and MIPS fluxes for all of the SINGS galaxies. Galaxies such as 
NGC 4725 or NGC 2841, with the lowest values of (vF^n^/ (yF^isQ^, require cooler dust (i.e., 
lower starlight intensities) than any of the SINGS-SCUBA galaxies. NCG 5408, with a very high 
value of (vF^-ji^/ (vF^iQQ^, requires hotter dust (higher starlight intensities). Submillimeter 
observations would be of great value to test the predicted submm fluxes from these extreme cases. 



^he 25 galaxies in Table 3 of foale et all (|2005h . plus NGC 1482 



2 



NGC 7552 is absent from Fig. 1 because of detector saturation at 71^m. 
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2.4. IRAS and ISO 

IRAS photometry is obtained from the latest SCANPI data for most SINGS sources, and from 
HiReS maps for the largest SINGS targets (optical diameters > 10'), except in some bands (e.g., 
no IRAS 12/im flux for NGC 0024). A 20% calibration uncertainty estimate was applied. In the 
case of NGC 5195 (=M51b), the proximity of NCG 5194 (=M51a) made flux determination more 
difficult, and 50% uncertainty estimates were adopted. 

Ths IRAS 25^m flux densities are not used, as they are superseded by the MIPS 24/um pho- 
tometry at essentially the same wavelength, but we employ the IRAS 12/im, 60/iin, and 100/iin 
photometry in our model-fitting. 

When available, we show ISOCAM 6.75/im and 15/xm photometry in the figures showing the 
spectral energy distributions (SEDs), but they are not used in the model- fitting. 



Physical Dust Models 



A direct estimate for the dust-to- hydrogen mass ratio Md us t/MH in the local Milky Way can 
be obtained from the difference between total interstellar abundances - taken to be equal to current 
estimates for solar abundances - and observed gas phase abundances in diffuse interstellar clouds, 
where the intensively studied diffuse cloud on the line-of-sight to COph is taken to be representative. 
Using current estimates of solar abundances, this mass inventory results in M^ nst /M}i ~ 0.0073 
(see Table [2]). 



We will apply (and ther eby also test) the d ust models developed by IWeingartner fc Drains 



(I2001al . hereafter WD01) and iLi & Draind (|200ll . hereafter LD01), and updated by DL07. These 
dust models consist of specified mixtures of carbonaceous grains and amorphous silicate grains, with 
the smallest carbonaceous grains having the physical properties of polycyclic aromatic hydrocarbon 
(PAH) particles; the size distributions are chosen to reproduce the observed wavelength-dependent 
extinction in the local Milky Way, in the Large Magellanic C l oud, and the Small Magellanic Cloud 
bar region ( Weingartner Sz Drainell2001al ). ILi &: Draind ( 200ll . l2002l ) showed that these grain models 
appeared to be consistent with the observed IR emission from diffuse clouds in the Milky Way and 
the Small Magellanic Cloud. 

The DL07 dust models used here adopt PAH optical properties that are slightly different from 
those adopted by LD01. As discussed by DL07, the changes are in part to better conform to recent 
theoretical estimates for ban d strengths, a nd in part to improve agreement with band profiles 
determined observationally by lSmith et al.l (|2007l ). 



The PAH emission features at 3.3/mi, 6.2/mi, 7.7/um, 8.6/um, 11.3/im, 12.0/um, 12.7/mi, 16.4/um, 
and 17 (im are prominent in many galaxy spectra, and the dust models considered here include 
various abundances of PAHs; in each model the PAHs have a wide range of sizes. Because PAH 
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Table 2. Dust Mass per H from Milky Way Abundances 



X (N 


x/N H ) Q (ppm) a (Nx/Nk 




M x , dust /M H 


C 


247 


0.57 


0.0013 


N 


85 


0.72 


0.0003 





490 


0.73 


0.0021 


Mg 


38 


0.08 


0.0008 


Al 


3" 


<o.i c 


0.0001 


Si 


32 


0.05 


0.0009 


Ca 


2" 


0.0002 d 


0.0001 


Fe 


29 


0.007 


0.0016 


Ni 


2 


0.004 


0.0001 


total 






0.0073 


a Jenkins ( 


2004) except as noted. 






b ( Nx/Nn 


)m from Grevesse & Sauval 


(1998) 





c assumed 



d Savage fc Sembach (1996) 



Table 3. Physical Dust Models 



JM 


Model 


?PAH 


M duBt /M^ b 


Av/Nh 






% 




mag cm 2 /H 


1 


MW3.1_00 


0.47 


0.0100 


5.3 x 10" 22 


2 


MW3.1.10 


1.12 


0.0100 


5.3 x 10" 22 


3 


MW3.1.20 


1.77 


0.0101 


5.3 x 10~ 22 


4 


MW3.1.30 


2.50 


0.0102 


5.3 x 10~ 22 


5 


MW3.1.40 


3.19 


0.0102 


5.3 x 10 -22 


6 


MW3.1_50 


3.90 


0.0103 


5.3 x 10~ 22 


7 


MW3.1_60 


4.58 


0.0104 


5.3 x 10~ 22 


8 


LMC2_00 


0.75 


0.00343 


1.2 x 10" 22 


9 


LMC2_05 


1.51 


0.00344 


1.2 x 10" 22 


10 


LMC2.10 


2.40 


0.00359 


1.2 x 10" 22 


11 


SMCbar 


0.010 


0.00206 


6.2 x 10" 23 


a M H 


= M(HI + H 


2) 







6 Mdust/M gas = [Mdust/M H ]/1.36 
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emission is mainly the result of single-photon heating, the PAH emission spectrum depends on the 
size (heat capacity) of the PAH. The prominent 6.2/iin, 7. 7//m, and 8.6//m em ission features are 



primarily due to PAHs with < 10 3 C atoms [see Fig. 7 of iDraine fe Lil (J2007J)]. Following DL07, 
we characterize the PAH abundance by a "PAH index" qpAU, defined to be the percentage of the 
total grain mass contributed by PAHs containing less than 10 3 C atoms: 

M(PAHs with N c < 10 3 ) 
4PAH = , (4) 

-Mdust 

where Md us t includes the mass of all silicate and carbonaceous dust, including the PAHs. 

WD01 put forward seven different size distributions for dust in the diffuse ISM of the Milky 
Way. The models are all consistent with the average interstellar extinction law, but have different 
PAH abundances, ranging from (?pah = 0.47% to 4.6%. Some properties for these seven models 
are given in Table [3j where the dust masses are calculated assuming densities of 3.5 gem -3 for 
silicates, 2.24 gem -3 for carbonaceous grains (WD01). The WD01 Milky Way dust models have 
-^dustA^H ~ 0.010, about 40% larger than the value 0.0073 estimated from observed depletions 
(see Table [2]), assuming solar abundances of C, Mg, Si, and Fe in the local interstellar medium. 
It is possible that current estimates of solar abundances of C, Mg, Si, and Fe may be lower than 
the actual elemental abundances in the interstellar medium today; a ~25% increase in the solar 
abundances would bring M<-j ust /MH estimated from depletions into accord with the WD01 model. 
Alternately, it is possible that the WD01 model overestimates the mass of interstellar dust by a 
factor 1.4 or so. 

The PAHs constitute the small size end of the carbonaceous grain population; the varying PAH 
abundances are the result of varying the grain size distribution, while holding the overall extinction 
constant. The model with the largest PAH abundance appears to reproduce the observed IR- 
submm emission from dust in the local interstellar medium (LD01, DL07). Because it is likely that 
the PAH abundance varies from galaxy-to-galaxy, here we take linear combinations of these models 
to produce a sequence of 43 MW models with ^pah = 0.4% to 4.6% in steps of 0.1%. 

Similarly, we take linear combinations of the 3 WD01 models (as modified by DL07) for dust 
in the LMC2 region (with gpAH = 0.75, 1.49, and 2.37%) to produce a sequence of 18 LMC models 
with gpAH = 0.7% to 2.4%. 

The dust models include amorphous silicate grains. We assume optically-thin dust when 
modeling the global SEDs, and therefore the models presented here do not include silicate absorption 
at 9.8/im. Such absorption is in f act seen in observ ed spectra of the nuclei of a few of the galaxies 



in this sample - e.g., NGC 3198 (jSmith et al J 120071 ) - but we neglect dust absorption in modeling 
the global emission. 

The dust models employed here do not include H2O and other ices. H2O ice is known t o be 



present in dark clouds in the Milky Way, in regions where Ay > 3.3 mag (jWhittet et al.lll988l ). In 
these regions, ices will modify the optical properties of the dust. However, it appears that such 
dark regions contain at most a modest fraction of the total dust in a galaxy; furthermore, the weak 
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starlight in such regions implies that the dust grains located there radiate only a very small fraction 
of the total IR power. As will be seen below, an ice-free dust model appears to be satisfactory for 
modeling the IR emission from normal galaxies. 



4. Constraining the Dust Model and Starlight Intensities 



In addition to choosing a physical dust mixture (relative numbers of carbonaceous grains and 
silicate grains of different sizes, including the fraction (?pah of the dust mass that is in PAHs), 
we must specify the intensity of the radiation that is heating the dust grains. Almost all of the 
heating of dust grains in a normal galaxy is from absorption of hv < 13.6 eV photons. In a normal 
star-forming galaxy a significant (but sub-dominant) fraction of the stellar luminosity may be at 
hv > 13.6 eV, but most of these photons end up photoionizing hydrogen or helium rather than 
being absorbed by dust. The hv < 13.6 eV recombination radiation - especially Lyman a - may 
be largely absorbed by dust. However, the total power in ionizing photons is a minor fraction of 
the total stellar luminosity in a galaxy with a more-or-less steady rate of star formation. Only in 
galaxies that have undergone an extreme burst of star formation within the last ~ 10 7 yrs will the 
dust heatin g be dominated by L yman a and Lyman continuum radiation. Only in rare cases, e.g., 
NGC 1377 (IRoussel et al.ll2006l ). where the massive stars formed in a starburst are surrounded by 
dense, dusty H II regions, will the dust grains directly absorb a major fraction of the hv > 13.6 eV 
radiation. 

The IR emission from dust is relatively insensitive to the detailed spectrum of the hv < 13.6 eV 
photons, and the DL07 models simply adopt the spectrum of the local interstellar radiation field 
as a reasonable representation for the average spectrum of the interstellar radiation in a normal 
spiral galaxy. The specific energy density of starlight is therefore taken to be 



u,, 



U X n (MMP83) 



(5) 



where lt ( / MMP83 ) j s ^he specific energy density estimated by Mathis et al. ( 1983 ) (hereafter MMP83) 



for the local Galactic interstellar radiation field, and U is a dimensionless scale factor. The intensity 
of interstellar UV radiation is often characterized by Go, the ratio of the 6-13.6 eV energy densit y 
relative to the value 5.288 x 10 -14 ergs cm -3 for the radiation field estimated by iHabingl ( 19681 ) . 
The MMP83 ISRF has u(0 - 13.6 eV starlight) = 8.65 x 1(T 13 ergscm~ 3 , and u(6 - 13.6 eV) = 
6.01 x xl(T 14 ergscm- 3 . Thus U = 0.88G . 

Dust at different locations in the galaxy will be exposed to different radiation intensities. A 
physically correct model would specify the locations o f the stellar source s and the absorbing dust, 
and would solve the radiative transfer problem [e.g., iTuffs et al.l (|2004l )]. but this would require 



many uncertain assumptions about the relative distributions of stars and dust, as well as massive 
numerical computations. 



Following DL07, we suppose that a large fraction of the dust in a galaxy is located in the 
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diffuse interstellar medium, exposed to starlight with an intensity that does not vary much within 
the galaxy. This "diffuse ISM" dust is idealized as being heated by a single radiation field, with 
intensity factor U = U m i n . In addition, DL07 suppose that a small fraction of the dust is located 
in regions where the radiation field is more intense, with a wide range of intensities ranging from 
U = U m - m to U = [/max ^> C/min- This would include dust in photodissociation regions, so we refer 
to this as the "PDR" component, and we use the power-law distribution function used by DL07. 
Although this "PDR" component contains only a small fraction of the total dust mass, in some 
galaxies it contributes a substantial fraction of the total power radiated by the dust. The models 
do not include a component corresponding to cold dust in dark clouds; as discussed in H9.1\ such 
cold dust appears to account for a minor fraction of the dust mass, and a very small fraction of the 
total dust luminosity. Thus, following DL07, we take the dust mass dM^ ast exposed to radiation 
intensities in [U, U + dU] to be a combination of a 6 function and a power-law, 

^§f^ = (1 - 7 )M dust 5{U - U min ) + 7 M dust J" ~ ^ U' a for U min < U < U max , a + 1 (6) 

where 5 is the Dirac 5 function. This distribution function has five parameters: the total dust 
mass Md U st, the fraction 7 of the dust mass that is associated with the power-law part of the 
starlight intensity distribution, and 3 parameters (t/ m in> ^max, ot) characterizing the distribution 
of starlight intensities in the high intensity regions. The 5-function term represents the fraction 
(1 — 7) of the dust that is heated by a general diffuse interstellar radiation field with ap proximately 



unifor m intensity ~ U m m- The power-law part of the distribution function, introduced by I Dale et al 



(|200ll ) and used by iLi &; Draind (|2002l ) , represents the fraction 7 of the dust mass that is close 
to OB associations and is exposed to radiation with intensity U > U m i n ; this includes the dust 
in "photodissociation regions" (PDRs), where radiation intensities can be orders of magnitude 
higher than the "average" dust grain is exposed to. Equation (|6j) is a very simplistic description 
of the distribution of dust grain environments; it neglects the possibility that the dust properties 
(composition, size distribution) might be variable and correlated with U, and also neglects regional 
variations in the spectrum of the starlight irradiating the dust. Despite its simplicity, we will see 
here that eq. © appears to capture the essential elements of the distribution of starlight intensities 
found in galaxies, yielding models with SEDs that approximate what is observed. 

Very small particles cool almost completely between photon absorptions, and therefore the 
shape of the emission spectrum from these grains is almost independent of U. The larger grains, 
that account for more than half of the total absorption of starlight in the Milky Way, do not 
cool significantly between photon absorptions, and their temperature depends on the intensity U. 
The observed FIR emission spectrum is sensitive to the temperature of these larger grains, and 
therefore allows us to constrain the distribution function for U. While Spitzer Space Telescope 
provides unprecedented sensitivity to FIR emission, the spectral coverage is somewhat sparse, with 
a factor of 3 gap in wavelength between IRAC 7.9^m and the MIPS 24^m bands, and a factor of 
2.9 gap between the MIPS 24/mi and 71//m bands. The IRS "Peak-Up" mode observations can be 
used for imaging at 16/um and 26/im, but these observations are not generally available. The large 
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gaps in wavelength coverage leave considerable uncertainty in the spectral energy distribution at 
wavelengths between the bandslf] Furthermore, the longest wavelength band, at 160/im, falls on the 
short wavelength side of the emission peak for dust temperatures < (l/6)/ic/(160 fim) = 15 K. As 
a result, Spitzer photometry does not strongly constrain the abundance of dust with temperatures 
Td < 15 K. To determine the amount of cool dust, global submillimeter observations are required. 
Fortunately, such data exist for 17 of the 65 galaxies for which we have complete Spitzer photometry 
(or near-complete in the case of NGC 7552), and this subsample will be analyzed first, in $5j 

We select a model for the IR emission from a galaxy as follows. For each size and composition 
of dust, we precompute and tabulate the temperature distribution function for the dust grain 
illuminated by the radiation field of eq. ([5]) for selected values of U (with 0. 1 < U < 10 7 ), usin g 



heat capacities, absorption cross sections, and numerical methods described by lDraine fe Lil (|2001 



LD01, and DL07. Neutral and ionized PAHs are assumed to have different absorption cross sections. 
The ionized fraction <fii(Nc) (a function of the number Nq of carbon atoms in the PAH) will 
depend primarily on the ratio U\/T P ~]n P , where n P and T e are the electron density and temperature 



( Bakes fc Tielens 1994 : Weingartner fe Drainel 



estimated by LD01 for the local diffuse ISIV^ 



2001b! ) . Here we take the ionized fraction <f>i(Nc) 



(we are therefore implicitly assuming that other 
galaxies have the same mix of Uy/Tl/n e as the Milky Way). Our model- fitting procedure assumes 
that all galaxies have the same shape for the PAH size distribution, the same PAH ionized fraction 
< j ) i{Nc)i and the same spectral shape of the illuminating starlight. As a result, our model fits have 
fixed ratios of PAH emission band strengths. Observations do reveal variations in PAH band ratios 
from one galaxy to the next - for example, L(7.7/mi)/-L(11.3/ii n) varies from ~4 -5 to <0.7 as we 
progress from H II-dominated nuclei to AGN-dominated nuclei (ISmith et alJl2007l ). and low values 
of L(7.7/im)/L(11.3//m) are also seen in some AGN- free galactic and extragalactic regions (e.g., 
Reach et al.ll2000l ; iHony et al.ll200ll : iPagani et al.lll999l ) - but most spiral galaxies have global PAH 
feature ratios that are near- average and characteristic of star-forming regions. 

With these temperature distribution functions for the dust, and the dust absorption cross 
section C a b s (A), we can calculate the time-averaged IR emission for a given grain type and size for 
each of the discrete values of U considered; we then sum over the grain types and sizes for a dust 
model j'm (see Table [3} to find the power radiated per unit frequency per unit mass of dust mixture 
ju exposed to starlight intensity U: 



1 



M(j M 



E 

k 



da 



dn k 
da 



J A I 



dr dP(M) 4 ^ 

al 



3 In principle, the IRS instrument can provide 6-36^m photometry, and the MIPS SED mode covers 55-95/im, 
but observing time considerations preclude using these observing modes for large areas. 

4 Following LD01, we assume 43% of the dust to be in the "cold neutral medium" (CNM; n e /U ~ 0.045 cm -3 , 
T « 100 K), 43% to be in the "warm neutral medium" (WNM; n e /U « 0.04cm" 3 , T w 6000 K), and 14% to be in 
the "warm ionized medium" (WIM; n e /U w 0.1 cm~ 3 , T w 8000 K). The resulting 4>i{N c ) is shown in Fig. 7 of LD01 
and Fig. 8 of DL07. 
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M(J M ) s eMIt) • < 8 > 

where the sum runs over the different grain types k (amorphous silicate, graphite, PAH , PAH + ) 
in the DL07 dust model, with p& the density of material k, and nfc(a) being the number of dust 
grains of type k per H nucleon, with radii < a. 

We also calculate the specific power per unit dust mass 

Pu(jM, Umin, U max , a) = \ _ > / p®(jM, U) U~ a dU (9) 

^min ^max JU min 

for a power-law distribution of starlight intensities, dM^ ust /dU oc U~ a . For each dust model jm 
we consider a range of values of U m \ n and U max . 

At wavelengths A < 5 fj,m, the emission from a galaxy is dominated by starlight. For A > 3.6 /im 
we neglect reddening, and appr oximate the s t ellar contribution by a dilute blackbody with color 



temperature T* = 5000 K, which ISmith et al.l (120071 ) found to provide a suitable approximation to 



the stellar continuum at A > 5/im. The model emission spectrum is then 

^.rnodcl = tt±B u (T±) + ^T^f (1 - 7)p1 0) {jM, U m in) + JPis(jM, ^min, ^max, Ol) , (10) 

where fi* is the solid angle subtended by the stars. For each (jM,U m - m ,U m & x ,a), we obtain 
Mdust, an d 7 by minimizing 

2 X - ^ [-^obs,6 — (-^,model)fe] 
X ^2^ ^2 +a 2 - (11) 

I) obs,b 1 model, b 

where the sum is over observed bands b, (-Fi/, m odel)& is t ne model spectrum convolved with the 
response function for band b, and cr b Si b is the observational uncertainty in the observed flux density 
-^obs.fi for band b. The ad-hoc insertion of the 0" modcl b term in the denominator in eq. (fTT|) is to 
allow for the approximate nature of our models and is motivated by the very high photometric 
accuracy achieved for the MIPS 24/im photometry, with <r b s ,24^m/-^obs,24^m = 0.040 achieved for 
many of the SINGS galaxies. We do not want the fitting procedure to give extreme weight to one or 
two bands where F b s ,& may have been measured with very high precision (i.e., small 0" o b s ,&/-fobs,&) 
because we do not expect our models to fit perfectly in the limit of perfect observations. Our 
subjective sense is that the present models should be regarded as satisfactory if they are able to 
reproduce the actual fluxes (i.e., perfect data, with cr bs,f> = 0) to within, say, ±10%. We therefore 
arbitrarily set <r mo dei,6 = Q-l(Fv,modei)b f° r each band b. It is important to recognize that eq. (fTTj) 
differs from the usual x 2 used for testing models, and the numerical values of x 2 found here are 
not rigorous measures of the goodness-of-fit. Nevertheless, these x 2 values do appear to provide a 
practical guide to distinguish between good fits and poor fits. 

The model components must be nonnegative: we require £1* > 0, M^ ust > 0, and < 7 < 1. 
We select the dust model jju, arid intensity model (7, U m { n , U max , a) giving the best fit (smallest 
X 2 ) subject to these constraints. 
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For the MW dust models, there are effectively 7 adjustable parameters: SI*, Md ust , (/pah (via 
3m) i 7, C^min, C^max, and a. However, we will fix one of these parameters from the outset, by setting 
a = 2. This choice is based on experimentation that indicated that the quality of the fit is not 
very sensitive to the precise value of a, with a « 2 providing satisfactory fits for diverse galaxies. 
Therefore we start with 6 parameters to be adjusted. For the galaxies with SCUBA data, we will 
examine the effects of varying U max , but we will see in §5.4.21 that setting £/ max = 10 6 works well; 
therefore, the number of adjustable parameters is, effectively, 5: $7*, Md us t> <7pah, U m \ n , and 7. 
If the number of observed bands is N b , the number of degrees of freedom = (iV& — 5) where Nj, 
ranges from N b = 9 (e.g., NGC 0024, with 4 IRAC bands, 3 MIPS bands, and IRAS 60, 100/xm) to 
N b = 12 (e.g., NGC 3190 with 4 IRAC bands, 3 MIPS bands, IRS 16//m, IRAS 12,60,100/im, and 
SCUBA 850/iin). To assess the goodness of fit, we calculate \ 2 P er degree of freedom: 

x ^i&i ■ (12) 



5. SCUBA-Constrained Dust Models 

5.1. Spectral Energy Distributions 

The distribution of grain temperatures in a dust cloud or a galaxy is never determined 
precisely. At wavelength A the thermal emission per unit dust mass oc [ e hc / XkT <i _ 1]~ 1 . For 
A < 360/im(20 K/Trf), small errors in the adopted lead to large errors in the dust mass inferred 
from the observed F v . For example, varying the adopted from 22 to 18K produces a 150% 
increase in the mass of dust required to produce a given 160/mi flux, but only a 35% increase in 
the mass required to produce a given 850/Lim flux. Observations at A > 450 ^m provide dust masses 
with the least sensitivity to dust temperature (for dust temperatures > 15 K). Therefore, the 
subset of 17 SINGS-SCUBA galaxies will be those for which we can best determine the dust mass. 

For each galaxy Table 0] lists the best-fit model for 10 3 < U ms , x < 10 7 , for fixed a = 2. If the 
best-fit model has either LMC or SMC dust models (see NGC 5195) then we also list the parameters 
for the best-fit model using MW dust. Finally, if the best-fit MW dust model has U mSLX / 10 6 we 
also list the parameters for the best-fit MW dust model with U m3iX = 10 6 . 



Table 4. 17 SINGS Galaxies with IRAC, MIPS, and SCUBA global fluxes 
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comp, q piH [U min ,U ma J 7 <U> lg(M d /M ) X ? E 

— MW 4.6% [2.5,lxl0 6 ] 1.44X10" 3 2.5 7.19 0.391 

— MW 3.8% [2.5,lxl0 7 ] 1.47xl0- 3 2.6 7.19 0.38- 



comp, q PAH [U min ,U mo J 7 <U> lg(M d /M e ) x? 1 

— MW 2.0%[10.0,lxl0 6 ]5.24xl0- i! 15.5 7.25 0.46 

— MW 2.4% [8.0, lxlO 5 ] 9.28xl0- a 14.3 7.29 0.45- 
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Fig. 3. — (a) SED for the SAap galaxy NGC 3190. Rectangles show observed fluxes in the IRAC, MIPS, SCUBA 
bands, and the IRS 16/im band; diamonds show 2MASS 2.2/im, and IRAS 12, 60, and 100/wn, fluxes. Vertical extent 
of rectangles and diamonds corresponds to ±la range, and width corresponds to nominal width of band. Solid line 
shows best-fit model U max = 10 6 model, with IRAC, MIPS, IRS 16/xm, IRAS, and SCUBA data used to constrain the 
fit; N b = 12 for NGC 3190. Triangles show model convolved with the IRAC, MIPS, IRS 16^m, IRAS, and SCUBA 
bands. Dot-dashed curves show separate contributions of starlight and emission from dust heated by U = U m i n 
(labelled "diff"), and dust heated by U m i n < U < C/ ma x (labelled "pdr"). Long-dashed curve shows best-fit model 
when C/max is unconstrained - for this case the model spectrum is nearly indistinguishable for the simple reason that 
there is relatively little dust (7 ~ 0.0015), and therefore relatively little power, in the "pdr" component, (b) As for 
(a), but for the SBa galaxy NGC 2798, with 7V 6 = 11 ( global imaging in the IRS 16/im band is unavailable). 

The models are generally very successful in fitting the observed SEDs. The SED of NGC 3190, 
together with a model fit, appears in Figure [3^,. This is an example where the model (with t7 m j n = 
2.5, [/max = 10 7 , 7 = 1.5 x 10~ 3 ) provides an excellent fit to the observations, with Xr — 0.38. 
The inferred dust mass is M dust = 1O 7 - 19 M , and the mean starlight intensity is estimated to be 
(U) = 2.6. Only a small fraction 7 ~ 0.0015 of the dust is in regions with starlight intensities 
U > C/ mm . Also shown in Figure [3^, (solid line) is a fit with the upper cutoff fixed at t7 max = 10 6 . 
The model SED is nearly identical (the solid curve and dash-dot curves coincide in the plot) and 
the derived dust mass is unchanged. 

Figure [3b shows the SED of NGC 2798, a galaxy where the dust is considerably warmer than 
in NGC 3190. The best-fit dust model (U min = 8, U max = 10 5 , 7 = 0.09) provides a good fit 
(Xr — 0.45) to the observations, but the model with t7 max = 10 6 provides nearly as good a fit 
(Xr = 0.46), a derived dust mass differing by only 0.04 dex, and similar (/pah (2.0% vs. 2.5%). 

The model fits to NGC 2798 both have a noticeable "bump" near ~ 33 /xm, produced by 
the warm graphite grains in the DL07 model. The small but nonzero conductivity for E \\ c in 



gra phite^ causes the absorption cross section C a \, s to rise to a broad local maximum near ~ 33 ^m 




5 The c-axis is normal to the sheets of hexagonally-organized C atoms. 
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log(A/ dust ) b 


log(Z, d ust) b 9pah c 


{U) d dust [/ min c (7 max ° 7 f 


/([/> 10 2 ) g 




type Mpc 


M Q 


Me 


M Q 


Le % 


type % 
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ngc5713 


Sbc 


26.60 
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10.02 


7.95 


10.72 


2.2 


6.63 


LMC 


4.0 
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26.8 


0.41 
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MW 
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0.54 












8.06 
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4.3 


3.20 


MW 


3.0 


10 6 


0.57 


4.9 


0.55 


ngc7552 


Sab 


22.30 


9.68 




7.84 


11.07 


1.5 


12.43 


MW 


8.0 


10 7 


4.24 


31.5 


0.50 












7.84 


11.07 


2.2 


12.55 


MW 


8.0 


10 6 


5.30 


31.1 


0.63 



a Z), 21 cm flux, and CO 1-0 flux from iKennicutt et al. I (|2003h unless otherwise noted. 
b Estimated from dust model. 

Traction of dust mass contributed by PAHs with Nc < W 3 C atoms. 
d Dust-weighted mean starlight intensity. 

c Umin, Uma, x = lower, upper cutoff for starlight intensity scale factor U. 
'Fraction of dust mass in regions with U > U m iu- 
g fraction of dust luminosity from regions with U > 10 2 
h X 2 /{Nb — 5), where Nt = number of bands used for fitting. 



'21 cm flux from lWalterl (|2005h . 
j CO 1-0 flux from lSheth et al.1 (|2005h . 



Freedman et all l|20011) . 
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( Draine &: Led Il984l . DL07). When exposed to starlight with U > 10 , graphite grains produce a 
broad emission bump peaking in the 30-35 /im region, depending on U. Unfortunately, the MIPS 
24/im and 71/mi bands straddle this feature, and do not constrain it. Implications of the model 
33/im feature are discussed below in £ 19.31 
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Fig. 4. — Histograms of Md us t/A / /(H I + H2) for the 12 galaxies in the sample for which both H I 21 cm and CO 1-0 



have been measured, using Xco =4x10 cm ( K km s 



to estimate M(H2) 



5.2. Dust-to-Gas Ratio 

If the interstellar abundances of carbon and all heavier elements were proportional to the gas- 
phase oxygen abundance, and if the same fraction of the major condensible elements (C, O, Mg, 
Si, Fe) were in solid form as in the Milky Way, then all galaxies would be expected to conform to 

M u (0/H) MW 

where (0/H)mw is the oxygen abundance in the local Milky Way, and the factor 0.010 is from the 
MW dust models in Table [3 

Table H] gives the estimated dust mass Md us t for the 17 galaxies. For many of the galaxies 
there are measurements of 21 cm emission and CO emission from which the masses of H I and H2 
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can be estimated (see Tabled]). The usual assumption is made that the H2 mass is proportional to 
the CO 1-0 luminosity: 

M(H 2 ) = 1.57 x 1O 4 M ( 9n X ™ _V-jL.V7 5co \ . (14) 

V ' ° V4 x 10 20 cm- 2 (Kkm s" 1 )- 1 ; V M P c y VJykms- 1 / 

where Xqq is the ratio of H 2 column density to antenna temperature integrated over the CO 1-0 
line, and Sco = M-o f F v (l — 0)du, where J F v (l — Q)dv is the CO 1-0 line flux. Use of the 
so-called CO 1-0 luminosity to estimate H2 masses remains suspect, and the value to use for the 
proportionality constant Xco is uncertain , with estimates for the Milky Way rangi ng from 3.0 x 



10 20 cm- 2 (Kkm s" 1 )" 1 (jYoung & Scovillelll99lh to 1.56 x 10 20 cm~ 2 ( K km s" 1 )" 1 (IHunter et al 
1997). He r e we use the value Xqo = 4 x 10 20 cm~ 2 ( K km s -1 ) -1 , the value recommended by 



Blitz et al.l (|2006l ) for giant molecular clouds in Local Group galaxies. We discuss the value of 
Xqo i n §9.41 where we argue that dust mass estimates for the SINGS galaxies favor Xqo = 
4 x 10 20 cm- 2 (Kkm s" 1 )" 1 . 

Figure H] shows the dust/hydrogen mass ratio for the 12 SINGS-SCUBA galaxies for which 
the total (atomic and molecular) gas masses are known. The derived dust-to-gas ratios are tightly 
clustered, with median M, dnst /M h ~ 0.005. The lowest dust-to-gas ratio is for NGC 2798, with 
M dnst /M H » 0.0023; the highest is for NGC 2976, with M dust /M H » 0.0089. 

The dust to gas mass ratios in Figure |4] are broadly comparable to the value ~ 0.007 estimated 
from gas phase abundances in the local Milky Way (Table [2]) or 0.010 estimated from models for 
the observed Milky Way extinction (Table EJ). The observed range 0.003 < M dust /M H < 0.01 is 
consistent with galactic metallicities between ~0.3x and ~lx solar, if a similar fraction of the 
heavy elements is in the form of dust as in the Milky Way. 

Based on Figure HI we conclude that the model-fitting procedure provides dust mass estimates 
for this sample of galaxies that are in reasonable agreement with the dust masses expected if 
depletions are similar to depletions in the Milky Way. Furthermore, the model appears to be able 
to reproduce the observed emission spectrum (see, eg., Fig. [3|). The metallicity dependence of the 
dust-to-gas ratio will be examined in §7.11 



5.3. Role of LMC and SMC Dust Models 

We now ask whether it is necessary to use the LMC and SMC dust models in the determination 
of dust masses - do the dust models developed for the LMC and SMC lead to improved fits to the 
observed spectra and, if so, how different are the inferred dust masses and PAH abundances? For 
each galaxy in the SINGS-SCUBA subset, we searched for the model that minimized % 2 ; for 16 of 
the 17 galaxies, a MW dust model was preferred, for NGC 5713 an LMC dust model was favored, 
but in no case did the SMC bar dust model give the best fit. For NGC 5713, where an LMC dust 
model gave the smallest Xri we a ^ so §i ye i n Tabled] the parameters for the best-fit MW model; the 
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derived Md us t differs from that obtained using the LMC model by only 0.01 dex, and (/pah = 3.1%, 
vs. 2.2% for the LMC model; when the MW model fit is constrained to have C/ max = 10 6 , the 
derived (/pah = 2.8%. Therefore, at least for the 17 galaxies with SCUBA photometry, the model 
fitting can be limited to MW-type dust models with minimal impact on the inferred dust masses 
or PAH index. 

In ^6] below we demonstrate that we can obtain good estimates for Md us t and (/pah even for 
galaxies lacking submm photometry, by limiting the search to models with C/ m in > 0.7. Here we 
compare estimates for (/pah and Md U st obtained using LMC and SMC dust models with those 
obtained using MW dust models. 



5.3.1. PAH Abundance (/pah 

Let <7pah(MW) be the best-fit estimates for (/pah when fitting with only MW dust models, and 
let qpah(^MC) be the estimate obtained using only LMC+SMC dust models. As discussed below 
(see §7.4p we exclude 4 galaxies for which (/pah cannot be reliably estimated. Because the LMC 
and SMC dust models are limited to (/pah < 2.4%, we further limit consideration to galaxies for 
which gp AH (MW) < 2.4%. The MW models have g PAH (MW) > 0.4%. Therefore for the 7 galaxies 
where the best-fit MW model has (/pah = 0.4%, we can conclude only that the actual (/pah < 0.4%. 

As seen in Figure [U the best-fit MW models and the best-fit LMC+SMC models generally 
have similar values of (/pah _ 16/19 galaxies have |</pah(MW) — (/pah(^MC)| < 0.4%. There is 
one extreme outlier: the SO galaxy NGC 5866. NGC 5866 has relatively little dust - the flux in 
IRAC band 3 is consistent with starlight alone, and there is only a modest excess due to dust 
in IRAC band 4. As a result, (/pah is not very securely determined for NGC 5866. When the 
fitting is limited to LMC and SMC dust models, the best fit is produced by the SMC model, with 
<7pah = 0.1%. However, this fit, with Xr = 1-92, is poor compared to Xr = 0.97 for the MW dust 
model, so the MW fit would be s trongly-favored over the SMC model. Spectroscopy of the center 



of NGC 5866 (jSmith et al.l 120071 ) shows strong PAH emission features, consistent with the value 
9pah(MW) rj 2.0%. 

NGC 6822 is another case of a large discrepancy between (/pah = 0.7% obtained from MW 
models versus the value (?pah = 0.1% obtained when fitting is limited to LMC, SMC dust models. 
Once again, the MW dust model produces a significantly better fit (x% = 1.13) than the SMC 
model (Xr = 2.04). 

For Tololo 89 (=NGC 5398), on the other hand, (/pah < 0.4% is obtained using the MW dust 
models, but the best-fit LMC, SMC model is an LMC dust model with (/pah = 1-1%. In this case 
the fitting is confused by the apparent lack of dust emission in IRAC band 3, despite a strong dust 
excess in IRAC band 4. The MW fit has xl = 2.50, and the LMC fit has xl = 2.47 - an indication 
that there is something wrong with either the data or the dust model. 



23 




Fig. 5. — (a) qpAH determined by fitting LMC or SMC models vs. (/pah determined by fitting MW models, for 
galaxies with (/pah(MW) < 2.4%. The estimates of (/pah are in good agreement: 16/19 galaxies have |(/pah(MW) — 
<?pah(:eMC)| < 0.4% (region bounded by broken lines), (b) Ratio of dust masses for best-fit LMC, SMC model and 
best-fit MW dust model, for galaxies with (/pah(MW) < 2.4%. The dust masses are in good agreement: the median 
ratio is 0.836, and 15/21 galaxies are within ±0.20 of the median. 

Since 16 of the 19 galaxies have (/pah(MW) — (/pah(^MC)| < 0.4%, we conclude that the 
estimates for (/pah are generally not sensitive to the details of the trial dust models, and restricting 
the fitting to the MW dust models will give reasonable results for (/pah- 

5.3.2. Dust Mass 

FigureEb shows, for the 21 galaxies with </pah(MW) < 2.4%, the ratio of dust masses obtained 
using the two different dust models. The median ratio is 0.84, with 15/21 falling within the interval 
M dust (LMC, SMC)/M dust (MW) = 0.85 ± 0.20. It is clear that the inferred dust mass is relatively 
insensitive to the details of the model dust mixture used in the modeling. 

Given that we obtain similar values for the dust mass M dust and the PAH abundance (/pah 
using either the MW or the LMC+SMC dust models, we will henceforth limit the modeling to the 
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MW dust models. This has the advantage of providing a set of models with qpAU ranging from 
0.4% to 4.6% (see Table [3j) allowing systematic comparisons of the level of PAH emission among 
different galaxies. The resulting estimates for dust masses, starlight intensity distributions should 
be relatively robust, and even the PAH abundances appear to also reliable to within Aijpah = 0.004. 
From this point on we will limit the model-fitting to the MW dust models. 
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Fig. 6. — (a) Distribution of minimum starlight intensity scale factor Umin- (b) Distribution of maximum starlight 
intensity scale factor C/ max . The fits were limited to MW dust models. 



5.4. Radiation Field Parameters 

5.4-1- Lower Cutoff Umin 

Figure [6^ shows the distribution of the best-fit lower cutoff U m \ n for the SINGS-SCUBA galax- 
ies. In no case did the best-fit model have U m i n < 1, even though the model-fitting procedure 
considered values as small as 0.1. We therefore conclude that - at least in the SINGS-SCUBA 
sample - only a very small fraction of the overall FIR emission is produced by cool grains within 
dark clouds (grains inside dark clouds, shielded from external starlight, will have (7 < 1 and will 
be cool, unless embedded star formation is present). The 160^m and 850/um emission appears to 
be produced primarily by dust in diffuse regions with U > 1. Cool grains in dark clouds could, in 
principle, contain a substantial fraction of the dust mass while contributing only a minor fraction of 
the total IR luminosity. However, we will see below ( §9. 1|) that the derived mass of dust heated by 
starlight with U > 0.7 is generally close to the value expected based on the observed total (atomic 



25 



~m i 1 i 1 i 1 i 1 i 

Mean Starlight Intensity for 
17 Galaxies with IRAC + MIPS + SCUBA data 



I i I 



median (U)=4.32 



o 




CO 


CD 






w 


CO 


cn 


CO 




CD 






m 


co 




CO 


r- 


CO 






LO 




CO 


£^ 




lO 






c- 




a 


a 


a 








a 


E 


co 




CD 


r- 


CD 


in 






co 




CO 


N 


w 


CO 


cn 




co 


CJ> 


03 


co 


CO 


CD 


m 






■* 


r- 


N 


o 


■* 


CO 


«sf 


10 






N 


a 


s 


f 


a 


| C 


f 


1 


c 


C 






4 




6 


8 


10 







0.8 1 



20 



Fig. 7. — Distribution of (U) for 17 galaxies with SCUBA data (fits limited to MW dust with U„ 
adjustable 7 and t7 m in-) 



10 b , with 



and molecular) gas mass, indicating that the model- fitting procedure is not overlooking substantial 
amounts of cooler dust in dark regions. 



5.4.2. Upper Cutoff (7 max 

Figure [6)3 also shows that the model-fitting procedure does not strongly favor any particular 
value of £7 max . We have experimented with fixing C/ max , and we find that fixing (7 m ax = 10 5 results 
in only an increase in x 2 P er g a l ax Y of 0.47, while the number of degrees of freedom per galaxy is 
increased by 1. Fixing U max = 10 6 works almost as well - the mean increase in \ 2 P er galaxy is 
0.70 - but fixing C/ max at 10 4 or 10 7 results in larger increases in average x 2 P er galaxy (in both 
cases > 1) relative to the fits where C/ max is allowed to vary. 

Because of the insensitivity of the results to the precise value of U max , we will simply set 
t^max = 10 6 and reduce the number of adjustable parameters from 6 to 5: $7*, Md ust , qpAH, U m i n , 
and 7. Having C/ max as large as 10 6 is not surprising. The ultraviolet intensity in the Orion 
Bar photodi ssociation region, for example, is ~ 3 x 10 times stronger than the local starlight 
background ( Marconi et alll998 ; Allers et al. hood ), and the dust within compact and ultracompact 



H II regions, and in the surrounding PDRs, is expected to be heated by radiation fields as large 
as ~ 10 6 . We therefore expect the IR spectrum of a star-forming galaxy to have a noticeable 
contribution from dust heated by U ~ 10 6 . From this point onward, we will discuss only MW dust 
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models restricted to J7 max = 10 6 . 



5.4-3. Mean Intensity (U) and Fraction of Power from High- Intensity Regions 

Let Pq be the power absorbed per unit dust mass in a radiation field U = 1. For the intensity 
distribution ([6]) with a = 2, the dust luminosity Ldust is obtained by integrating 



dL dust = U P dM dnst = P M dust 
Thus the total dust luminosity is 



(1 _ l)u5{u _ Umin) + 1 u- 1 

U min ^max 



dU . (15) 



L dust = (U)P M dnst , (16) 
where the dust-weighted mean starlight intensity scale factor is 

7m([/ max /[/ mm ) 



(U) 



(i-i)u min + 



rr 1 - ij- 1 



(17) 



Figure shows the distribution of (U) for the 17 SINGS-SCUBA galaxies, fitted with MW dust 
models with U max = 10 6 . These galaxies are clearly not all alike: (U) varies by a factor of 6 over 
the sample, from (U) = 2.6 (for NGC 4569) to (U) = 16 (for Mrk 33), with median (U) = 4.3. 

The fraction of the dust luminosity radiated from regions with U > U c (for i7 m i n < U c < U max ) 

is 

fir .TT^TT\- iHUm^/Uc) 

J{^dust,U>U c —— -— — — — — — — . [IS) 

(1 - 7)(1 - C/ min /C/ m ax) + 7l n (^max/^min) 

Figure [8] shows the distribution of f(L dust ;U > U c ) for U c = 10 2 and 10 3 for the 17 galaxies with 
SCUBA data. The median galaxy in the SINGS-SCUBA sample has ~ 10% of the dust power 
originating in regions with U > 10 2 . In extreme cases, this fraction can be considerably larger, e.g., 
Mrk 33, with f(L dust ;U > 10 2 ) » 48%, and > 30% for NGC 2798 and NGC 7552. At the other 
extreme, the best-fit model for NGC 5866 (see £ }5im below) has 7 = and f(L dust ; U > 10 2 ) = 0. 
Values of U > 10 2 are expected to arise primarily in star-forming regions where dust is found 
near luminous stars. The values of f(L dus t, U > U c ) estimated from this fitting procedure are 
only approximate, but values of f(L dust ;U > 10 2 ) > 0.03 should be indicative of significant star 
formation rates. 



5.5. PAH Abundances 



Observatio ns with ISO show ed that normal star-forming galaxies routinely have strong 6-9/iin 
PAH emission (IDale et al.ll200ll ). This is confirmed for the present sample of 17 SINGS-SCUBA 
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Fig. 8. — Distribution of the fraction of the dust luminosity radiated by dust exposed to starlight intensities U > 10 2 
and U > 10 3 . MW dust models are used for the fits, with fixed (7 max = 10 6 , and adjustable 7, U m i n . Galaxies with 
/ < 0.01 have been grouped in a bin at 0.01. 



galaxies. Figure O shows a histogram of ^pah for these galaxies. All of the galaxies show PAH 
emission, with median iJpah = 3.5%. None of these galaxies have (/pah < 1-3%. The distribution 
of PAH abundances over a total of 61 SINGS galaxies will be discussed in ^7.41 below. 



5.6. Discussion of Selected Galaxies 

5.6.1. Discrepant Cases: NGC 3521 and NGC 5866 

Most of the 17 galaxies in the SINGS-SCUBA sample have SEDs that are satisfactorily repro- 
duced by the adopted dust model (see Fig. [3]). If we consider only MW dust models, with fixed 
kmax = 10 6 , the 17 galaxies have median Xr = 0.55; 16/17 of the galaxies have Xr < 1- The 
SINGS-SCUBA galaxies with the worst fits (largest xl) are NGC 3521 ( X 2 = 1.37), and NGC 5866 
(Xr = 0.97). Here we consider why the model has difficulty fitting these two galaxies. 

In the case of NGC 3521, the best-fit model underpredicts the 7.9/im emission, but the more 
significant discrepancies are at long wavelengths, where the model is high at 71/im, low at 160/mi 
and high at 850/im. To illustrate the difficulty with reproducing these three fluxes, Figure [T0l also 
shows a single-temperature component with a v 2 opacity with the temperature and power chosen 
to fit the 71 and 160^m data. Note that even this single-temperature component is too "broad" - it 
exceeds the 850^m datum. If a range of dust temperatures is present, it will broaden the spectrum, 
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and worsen the fit. Either there is a problem with some of the photometric data, or the dust in 
this galaxy has an opacity that is steeper than v 2 . 

The second-worst fit is for NGC 5866, an SO galaxy with a substantial edge-on disk, and the 
lowest value of (vF v )i±/ (vL v )yi among the SINGS-SCUBA galaxies (see Fig. [[]). The best-fit model 
has 7 = - all of the dust is exposed to starlight with intensity U = U m { n = 5. The small value 
of 7 is presumably because the illumination of the dusty disk is dominated by the distributed old 
stellar population, with little contribution fr om star formation in t he disk. This is consistent with 
nondetection of Ha, with SFR< 0.1M Q yr _1 (jKennicutt et alj|2003h (altho ugh dust extinction may 
also be important). The nuclear spectrum of NGC5866 (jSmith et al.l 120071 ) includes ionic emission 
lines (e.g., [Ne II], [Si II], [S III]) but they are relatively weak, consistent with only low levels of star 
formation; the low ratio [S III] /[Si II] flux ratio also suggests a deficiency of massive stars. 



Roberts et al.1 ( 199ll ) estimated the dust mass to be 1.4 x 1O 6 M0 (for our adopted distance 
of 12.5 Mpc). Here we estimate the dust mass to be 4.5 x 10 6 M Q - larger by a factor three. Our 
estimate for M dm t resul ts in M dust /M H ~ 0.005, consistent with A Q = log 10 (O/H) + 12 = 8.43±0.18 
dMoustakas et al.ll2007h . ~55% of solar ((A )© = 8.69 from Table©. 
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Fig. 10. — (a) SED of the SABbc galaxy NGC 3521, with xl = 1-37 for the best-fit MW dust model with 
C/ max = 10 6 . The dotted curve shows emission from hypothetical T = 21.5 K dust with ft„ oc v , showing that dust 
with k v oc v 2 cannot reproduce the observed 160/^m emission without exceeding either the 71/^m flux (if T > 21.5 K) 
or the 850/^m flux (if T < 21.5 K). Curves labelled "diff" and "pdr" show contributions of dust heated by U — U m i n 
and U > Umin, respectively, for the MW model with U max = 10 6 . (b) SED for the SO galaxy NGC 5866. The best-fit 
MW dust model - with a single starlight intensity (7 = 5- has Xr = 0.97. The model has difficulty reproducing the 
sharpness of the peak in \F\ near lOO^im and overestimates the flux at 450/im. 



Once again, the difficulty is in fitting the long-wavelength data. In this case, a single- 
temperature T = 24.5 K component with a v 2 opacity can fit the 71, 160, and 850/iin photom- 
etry, but exceeds the 450//m datum. Once again, a distribution of temperatures would exacerbate 
this problem, indicating either problems with the photometry or requiring an opacity substantially 
steeper than u 2 between 160 and 450/im. 



5.6.2. Markarian 33 



The dwarf starburst galaxy Mrk 33 (see Fig. ITT]) has the highest value of (vF v )2±l ' {vF v )i\ among 
the SINGS-SCUBA galaxies (see Fig.QJ. Ha observations of Mrk 33 indicate a SFR« 1.5M yr _1 , 
giving it the third-highest SFR/L opt ratio in the SINGS sample (jKennicutt et al.ll2003l ). The best-fit 
^max = 10 6 model has a very high mean starlight intensity (U) ~ 16, and /(Aiust! U > 10 2 ) « 48%, 
the highest value in the sample. 



6. Estimating Dust Masses Without Submm Data 

Thus far we have been fitting dust models constrained by IRAC, MIPS, IRAS, and at least 
one submm flux measured by SCUBA. Unfortunately, submillimeter observations are unavailable 
at this time for most of the SINGS galaxies. How well can dust masses be estimated in the absence 
of submm constraints? 
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Fig. 11. — Observed fluxes and model emission spectra for galaxies with IRAC, MIPS, and IRAS data, 

and at least one submm flux measured by SCUBA. Some galaxies (e.g., NGC 3627) have also been measured in the 
IRS 16/im peakup filter. Rectangles and diamonds show observed fluxes; vertical extent is ±lcr. Curves show models 
fitted to IRAC, MIPS, SCUBA, and IRAS 12, 60, and 100/im data,; triangles are models convolved with bandpass 
used for observations. Solid curves are models with MW dust and Umax = 10 6 ; broken curves are used for other 
models that for some galaxies give smaller Xr ( see text). Dashed-dotted curves show separate contributions of stars, 
dust with U = LUn (labeled "diff") and dust with U > t/ mill (labeled "pdr"). See Fig. for NGC 2798 and 3190, 
and Fig.[T0]for NGC 3521 and NGC 5866. 



To address this question, we first repeat the dust modeling procedure used in £JS] for the 17 
SINGS-SCUBA galaxies, but with the SCUBA data omitted. The dust masses so obtained are 
shown in Figure [121 where they are plotted versus the dust masses estimated when SCUBA data 
are also used to constrain the models. 
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Fig. 1111 continued. 



If the dust masses obtained with the SCUBA fluxes included are regarded as our "gold stan- 
dard", it is clear that dropping the SCUBA data from the fitting procedure leads to substantial 
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M dust ( M ©) [SST + IRAS + SCUBA] 

Fig. 12. — Dust mass Md us t determined for 17 galaxies using IRAS and SST data only, vs. the masses derived 
using IRAS, SST, and SCUBA data combined. Data are fit by MW dust models with [/max = 10 6 but no restriction 
on J7min- Without SCUBA data, cool dust is not strongly constrained. Nevertheless, 11/17 of the galaxies fall within 
a factor 1.5 of the value obtained when SCUBA data are employed, and all 17 galaxies are within a factor of 2.2 . 

uncertainties in the model fitting. Five of the 17 galaxies have estimated dust masses that exceed 
the "gold standard" estimate by more than a factor of 1.5 (NGC 3521, NGC 3627, NGC 4536, 
NGC 4569, and NGC 7331). In five cases the dust mass estimate is lowered when the SCUBA data 
are not used, in one case by a factor of more than 1.5 (NGC 5195). 

Figure [13] shows estimated dust-to-gas mass ratios for the 17 galaxies with SCUBA data for 
which the gas mass is known. Figure [T3"r shows the dust /gas ratios estimated when the SCUBA 
fluxes are used in the fitting: the median dust /hydrogen mass ratio is 0.0053 for this sample, 
comparable to the value ~0.007 expected for solar-metallicity material with ~2/3 of the C, and 
most of the Mg, Si, and Fe incorporated into carbonaceous and silicate material (see Table [2]). 
Figure [13b shows the effect of fitting the SED without using the SCUBA data - for most of the 
galaxies, the estimated dust/hydrogen ratio rises modestly, with a median value of only 0.0077. 

Because these models reproduce the SED fairly well, and the IR power is mainly shortward 
of 200 /im, the total IR power is relatively unaffected when the SCUBA data are removed. What 
is changing is the estimated mean dust- weighted starlight intensity, with (U) often falling a factor 
~ 1.5 or more below the value obtained when the SCUBA data are employed. In the absence of 
submm measurements to constrain the mass of cool dust, there is a risk that the model- fitting 
procedure may invoke a large mass of cool dust, heated by weak starlight. To prevent this, when 
submm data are unavailable we will restrict the model-fitting to use radiation fields with U m i n > 0.7. 
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This of course runs the risk of underestimating the dust mass if there are galaxies without submm 
(e.g., SCUBA) observations which actually contain large amounts of cool dust heated by starlight 
with U < 0.5. However, we note that all the galaxies with SCUBA data had U m i n > 1 for the 
best-fit model (see Fig. [6^). 
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Fig. 13. — (a) Histogram of A/ dus t/M H obtained using IRAC+MIPS+SCUBA data for 12 galaxies for which 
H I 21cm and CO 1-0 fluxes have been measured, (b) Same as (a), but estimated without using SCUBA data, with 
(7 max = 10 6 but no restrictions on C/ m i n - 



Based on the above discussion, we adopt the following "restricted" fitting procedure when 
submm data are unavailable: 



1. Models limited to MW dust models only, with 0.4 < g PA H < 4.6%; 
2 77 — 10 6 - 

3. a = 2; 

4. 0.7 < U mhl < 25. 



This restricted fitting procedure (with U m \ n > 0.7) of course implies that the resulting fits will 
necessarily have (U) > 0.7. This is generally consistent with what we observe for the sample of 17 
galaxies with SCUBA data - as seen in Figure [71 these galaxies tend to have (U) > 2. However, 
the SINGS-SCUBA galaxy sample may be biased in favor of increased dust mass and associated 
star formation. If other galaxies have very low mean starlight intensities (U), the restricted fitting 
procedure adopted here would underestimate the dust mass (because it would overestimate {U}) 
To study the dust in such galaxies, submm observations are essential. 
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7. Dust Properties for 65 SINGS Galaxies 

7.1. Dust to Gas Mass Ratio 

For 48 galaxies in the SINGS sample we have complete IRAC and MIPS data, but no SCUBA 
fluxes. Dust masses have been estimated for these galaxies using the restricted fitting procedure 
described above. The SEDs for these galaxies, together with the best-fit models, are given in 
Appendix [A] (Fig. [24"|) . The results for the dust mass Md us t are given in Table [5j 

For 20 of these galaxies both H I 21cm and CO 2.6mm fluxes have been detected (or there is 
a strong upper limit on the CO flux), and we are able to estimate the total gas mass. For another 
24 galaxies H I 21cm has been measured, but the CO 2.6mm flux is unknown, so we have only a 
lower limit on the gas mass. For 4 galaxies (the E galaxies NGC 0855 and NGC 4125, and the SO 
galaxies NGC 1266 and NGC 1316) neither H I nor CO have been detected. 

Figure E] shows the distribution of M^ ust /Mfi for the 20 galaxies lacking SCUBA data for 
which Mh is known. The distribution looks similar to the distribution in Figure 113b . except for 
two galaxies - IC 2574 and NGC 4236 - with M dust /M H ~ 0.0005. Note, however, that many of 
the sample galaxies do not appear on the histogram because Mh is not known. 

The dust-to-gas mass ratio reflects the enrichment of the gas by heavy elements (C, O, Mg, Si, 
Fe) produced in stars, the production and growth of solid grains from these elements in supernovae, 
stellar winds, and in the interstellar medium, and destruction of grains by supernova blastwaves 
and other violent events in the interstellar medium. 

If the same fraction of the condensible elements (C, Mg, Si, Fe, ...) is in grains as in the local 
Milky Way, then Md ust /MH should conform to eq. (|13p . However, it is possible that in some other 
galaxies, the balance between grain formation and destruction may be such that most interstellar 
C, Mg, Si and Fe is in the gas rather than in grains. We therefore expect eq. (1131) to define an 
upper envelope to the dust-to-gas mass ratio. 



Table 5. 48 SINGS Galaxies with IRAC and MIPS global fluxes, but no SCUBA data 1 
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a All models assume MW dust, with f7 max = 10 . 

b D, 21 cm flux, and CO 1-0 flux from Kennicutt et al. (2003) unless otherwise noted. 

c M(H a ) calculated using X C o = 4 x 10 20 cm- 2 (Kkm s" 1 )" 1 

d 21 cm flux from Walter (2005) 

c CO 1-0 flux from Sheth et al. (2005) 

f Tosi et al. (2001) 

g D is weighted mean of Ford et al. (1996) and Aihar et al. (1997) 
h CO 1-0 flux from Baiaia et al. (1991) 
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Models have be en proposed to attempt to predict the dust-to-gas ratio in galaxies as the y 
evolve in me t allicit y (|Dweklll998l ; iLisenfeld Ferraralll998l ; lEdmundsll200ll ; iHirashita et alj|2002l ). 
James et al.l (|2002l ) compared models to dust masses estimated from IR and submm observations, 
and argued that the fraction of the metals condensed into dust grains appeared to be the same in 
dwarf and giant galaxies - i.e., dust-to-gas mass ratios consistent with eq. (fl~3j) - for 8.1 < Aq < 9.0. 



~i — i — i i i i | 1 1 — i — i — i i i i | 1 1 — i — i — i i i i 
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Fig. 14. — Dust to H mass ratio determined from IRAC + MIPS data only, using restricted fits with U m i n > 0.7, 
and taking X C o = 4x IO 20 cm^(Kkm s -1 ) -1 . 



Galliano et al.l (|2003l . 120051 ) estimated the dust masses for 4 dwarf galaxies, and found that 
NGC 1569 and H e 2-10 had M dust /M H a factor of 4-7 and 2-10, respectively, below eq. (fT3l) . 



Hunt et al.l (|2005l ) used IR and submm observations to estimate dust masses for a sample of 
seven low-met allicity blue c o mpact dwar f s. Fo r SBS 0335-052, with metallicity Ao = 7.29 ± 
0.01 (jlzotov &: Thuanl Il998l ). iHunt et al.l (|2005l ) found an extremely low dust-to- gas mass ratio 
Mw/Mh ~ 1Q~ 7 , far b elow the value ~ 2 x 10~ 4 given by eq. (fl~3j) . The resu lts oflGalliano et al 



( 2005 ) and Hunt et al. ( 2005 ) are clearly inconsistent with the conclusion of James et al.l ( 20021 ) 
that dwarf galaxies have the same dust-to-metals ratio as the Milky Way. 



Figure [T5"a.b shows how M^ ust /M}i correlates with oxygen abundance for the SINGS galaxy 



Table 5 — Continued 



galaxy morph 


D h 


log[M(HI)] b 


Iog[Af(Ha) 


] b ' c log(A/ dust )° 


l0g(Ldust)° 


qpAH p {U) q dust 


t/ min r 7 s 


f(u>w 2 y xl u 


type 


Mpc 


M 


M 


M 




% type 


% 


% 



1 21 cm flux for radius < 6' from ^oulanger & Viallcfond (1992) 

j CO 1-0 flux from Lerov et al. (2005) 

k ?pah very uncertain - see Q7A\ 

1 Karachentsev et al. (2002a) 

m Karachentsev et al. (2002b) 

" Cannon et al. (2005) 

° Estimated from dust model. 

p Fraction of dust mass contributed by PAHs with Nc < 10 C atoms. 

r Lower cutoff for starlight intensity scale factor U. 

B Dust-weighted mean starlight intensity scale factor. 

t Fraction of dust luminosity from regions with U > 10 2 

u X 2 /(Nb — 5), where Nb = number of bands used for fitting. 
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Fig. 15. — Ratio of A/dust to Mh vs. oxygen abundance (see text) for galaxies with and without SCUBA 
fluxes (circles and diamonds, respectively), showing 32 galaxies where gas mass is known, for two values of 
Xco/[cm~ 2 (Kkm s" 1 )" 1 ]: (a) 2 x 10 20 , and (b) 4 x 10 20 . Filled diamonds: M dust /M H for regions where IR 
emission is detected (see text). Solid line shows eq. (|13j l; broken lines show factor of 2 variations around the solid 
line. Histogram of normalized dust-to-gas mass ratio for Xco/[cm -2 ( K km s -1 )" 1 ] = 2 x 10 20 (c), and 4 x 10 20 (d). 



sample, where we use the "cha racteristic" oxyge n abun dance j§ from Moustakas et al. ( 2007 ). based 
on only strong nebular lines. iMoustakas et al.l (|2007l ) show that different "strong-line" methods 



6 The characteristic abundance is intended to be represe ntative of the luminosi ty-weighted ISM material. For 
normal spirals, it is taken to be the metallicity at OAR25 (jMoustakas et al.ll2007l ). The quoted uncertainties in 



oxygen abundances do not include systematic uncertainties in the calibrations. 
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Fig. 16. — Galaxies with and without SCUBA fluxes (circles and diamonds, respectively), including 28 galaxies 
with only lower limits on the gas mass. Solid line shows eq. (|13[) ; broken lines show factor of two variations around 
the solid line. Note that the upper bounds are all galaxies where Md ust has been determined, but either H I or CO 
has not been observed. Filled diamonds show Mj DSt /Mh for regions where IR emission is detected (see t|7.2jl . 



[e.g. jKobulnickv fc Kewlevi (12004r) vs. IPilvugin fc Thuanl (|2005l )] give differi ng estimates for (O/H), 
with t he IPilvugin &: Thuanl (|2005l ) values tending to be smaller than the iKobulnickv fc Kewlev 



2004 ) values by a factor of four, on ave rage. Here we adopt abu ndances based on the IPilvugin fc Thuan 



20051 ) method, and refer the reader to lMoustakas et al.l (|2007l ) for discussion of the various system- 
atic uncertainties. For the local MW we take the gas-phase Ap = 8.59 as determined by absorption 



line studies of low-density sightlines (jCartledge et al.l 1 2 0041 ) . Figures [T5a and b differ only in the 
assumed value of Xqo, used to estimate the mass of H2. For both trial values of Xco, galaxies 
with 8.1 < Aq < 8.5 appear to conform to eq. (fT3j) to within a factor ~2. 



The two lowest metallicity galaxies in Figure El IC 2574 and NGC 4236, both fall well below 
eq. (|13|) , seemingly showing these systems to have only a small fraction of their condensible elements 
in dust. However, we will see below that this is not necessarily the case. The points labelled i2574* 
and n4236* show estimates for the dust /gas ratio for the regions in these galaxies where IR emission 
is observed (see ^7.21 for further details). IC 2574 now falls right on eq. (fT3|) . NGC 4236 is low by a 
factor ~ 2.5, but the estimate for the mass of gas in the IR-emitting region is very uncertain (see 



The highest metallicity system in our sample is NGC 4594 (the Sombrero galaxy, M104). 
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NGC 4594 also has the highest dust-to-metals ratio; even for Xco = 4 x 10 20 cm~ 2 (Kkm s" 1 )" 1 , 
(M,j us t/MH)/(0/H) is a factor of ~4 above eq. (JT3J) . Although this galaxy harbors a m assive AGN, 
the dusty disk and ring are strongly detected at 24, 71, and 160/im (jBendo et al.ll2006al ). The dust- 
to-gas ratio estimated for the disk of this galaxy is quite high, Md us t/AfH = 0.05; the reason for this 
high value is not clear. Perhaps the dust properties in this system are significantly different from 
the Milky Way dust properties that have been assumed in the modeling. Alternatively, perhaps the 
molecular mass in the dusty ring has been underestimated. The dusty ring has only been observed 
at a few positions in CO; additional CO observations would be valuable. It is also possible that 
the ratio of CO luminosity to molecular mass in this galaxy may differ from that in spiral galaxies 
where the so-called "X factor" has been calibrated. Unfortunately, the dust mass estimate for 
NGC 4594 is not yet constrained by submm data - while NGC 4594 is strongly detected at 850/um, 
the flux is dominated by the A GN, and the 850//m flux from the dusty disk has not yet been 
determined (jBendo et al.ll2006al ). We have seen above (Fig. [T2"j) that dust mass estimates made 
with and without SCUBA data can differ by factors of ~2. Measurement of the submm continuum 
emission from the dust ring in NGC 4594 would be of great value. 



7.2. Dust-to-Gas Mass Ratios in Low-Metallicity Galaxies 



There are many galaxies for which we have only lower limits on the gas mass, because either 
HI or CO has not been observed. The resulting upper limits on M^^/M-^ are included in Figure 
TB"| which now includes systems with Ao as low as 7.54±0.35 (Hoi). For Ao > 8.1, upper limits on 
the global dust-to-gas ratio are generally consistent with eq. f)13[) . to within a factor ~2. However, 
every galaxy in Figure [TBI with Aq, < 8.1 has a global dust-to-gas ratio that falls below eq. (|13p . 
sometimes by a very large factor. IC 2574 is low by a factor of 5, and 4 galaxies (Ho II, DDO 053, 
NGC 5408, NGC 2915) are factors of 8 or more below eq. CE]). The most extreme case is NGC 2915, 
below eq. (]13p by a factor of 40. At face value, this might suggest that the balance between dust 
grain formation and destruction in low-metallicity galaxies such as NGC 2915 is much less favorable 
to grain formation and survival than in the Milky Way. However, we will argue here that (1) the 
dust-to-metals ratio in the central regions of these galaxies is in fact normal (i.e., consistent with 
eq. [13]) . and (2) the dust-to-metals ratio outside the central regions is not yet known. 



Walter et alJ (|2007l ) compare the stellar, IR, and H I morphologies of the dwarf irregular (Im) 



galaxies Ho II, Ho I, IC 2474, and DDO 053: in each of these galaxies, most of the H I is located 
outside the region where IR emission is detected. For example, only a fraction /(H I) ~ 0.19 of the 
H I in IC 2574 lies within the regions where 70 and 160/im emission is seen. Similarly, the fraction 
of the H I within the contours where dust emission is detected is /(H I) 0.2, 0.29, and 0.33 for 
Ho II, Ho I, and DDO 053. 

We have estimated the radii #tr containing the bulk of the infrared emission for four additional 
low-metallicity dwarfs : NGC 1705 (|Cannon et alihoOBah . NGC 2915, NGC 4236, and NGC 6822 
(jCannon et alJ l2006bl ). Our estimates for and for the fraction /(H I) of the H I located at 
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9 < 6>i R , are given in Tabled with / ranging from 0.028 (for NGC 2915) to 0.5 (for NGC 4236). 
The dust-to-gas mass ratios for the regions where IR emission is detected are shown in Figure [T6j 
Six of the eight galaxies in Table [6J now have upper limits consistent with (|13|) to within a factor 
two, with only two (DDO 053 and NGC 4236) falling low by a factor ~2.5. Given the uncertainties 
in estimation of both M^ ust and /(H I), it is entirely possible that > 50% of the interstellar Mg, 
Si, and Fe may be in solid grains, at least in systems with Ao > 7.5. 

What can we say about the dust/metals ratio in the extended H I envelopes? Both the metal 
abundances and the dust masses in these regions are uncertain. The metallicities are estimated 
using H II regions, and therefore apply only to regions with recent star formation. The outer H I 
envelopes lack bright H II regions, and therefore the metallicity there is unknown. It would not 
be surprising if it were subtantially lower than the "characteristic" metallicities used in Figure 
[TBI For the outer regions of some spiral galaxies, oxygen has been found to be underabundant, 
relative to the "characteristic" metallicity of the central regions, by factors of four or more (e.g., 



Ferguson et al.lll998l : iTiillmann et aU 12003 ). 



In the case of NGC 1705, UV absorption line studies show that the H I envelope has Ao = 
7 A3 ± 0.22 (jHeckman et al.ll200ll ): the oxygen abundanc e in the H I envelope i s a factor 6 below 
the oxygen abundance determined from the H II regions (jMoustakas et al.l 120071 ). For this oxygen 
abundance, Mdust/Mn for NGC 1705 is consistent with eq. (fHfl) . 

The dust content of the envelope is often also uncertain. The stellar densities in the outer 
regions of these galaxies are low, and the starlight will be dominated by light coming from the 
central regions. Table [6] lists (vL v )b, the total stellar luminosity in the B band. For each of these 
galaxies we give the estimated half-mass radius #o.5 for the HI, and we have estimated the starlight 
intensity scale factor U at the H I half-mass radius, 

1 {vL v ) B 



U, 



0.5 



(1.982 x 10- 13 ergscm- 3 )47Tc(L>6» . E 



(19) 



where 1.982 x 10~ 13 ergs cm -3 is {vu v ) b for the local interstellar radiation field ( Mathis et al.lll983l ). 
Equation (|19j) neglects internal extinction within the dwarf galaxy [illumination of the outer en- 
velope will suffer the same attenuation as radiation escaping the galaxy, from which {uL v )-q was 
obtained]. Only one galaxy in Table 0- the dwarf starburst NGC 1705 - has [/0.5 > 0.2. 

The a > 0.03/^m grains that normally dominate the 160/mi emission are expected to have 
temperatures T « 17U 1 / 6 K = 12(f7/0.1) 1 / 6 K, with peak up v at A ss hc/QkT » 180(0.2/J7) 1/6 /im. 
Therefore, for the seven galaxies in Table[6jwith [To. 5 < 0.2, even the MIPS 160/um band (with half- 
response points at 140 and 174//m) is relatively insensitive to whatever dust may be present at and 
beyond #0.5 except for a small fraction that may be close to local stellar sources. MIPS photometry 
for these galaxies is dominated by the dust in the luminous regions, and the SED fitting is therefore 
fitting the parameters for this dust. If there is dust present beyond the half-mass radius, submm 
observations will be required to detect it in all of the galaxies in Table [6] except NGC 1705. 



In the case of SBS 0335-052, where iHunt et al.l (|2005l ) report an extraordinarily low dust-to-gas 
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Table 6. Parameters for Selected Dwarf Galaxies 
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Becker et al 



( 198£ 



b Estimated from Rotsl ( 1980[ ) map, assuming FWHM=10' gaussian beam. 



Walter et alj (|2007r i 



d Estimated from H I data in lMeureri (|1997T ) 
c Estimated from H I profile in lCannon et al.1 1 2006bh 
Meurer et al. |l998h 



j Cannon et al. (2006a 
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mass ratio Mj ust /MH = 1 x 10~ 7 , a factor 5000 below eq. (|13p . we note that this galaxy has not been 
detected at wavelengths between 100/mi and 1 cm. The dust models are poorly constrained, and 
it is possible that they may have significantly underestimated the mass of cool dust. In addition, 
much of the HI mass in SBS0335-052 may be in an extended envelope, where the metallicity is 
unknown and where the dust temperatures would be very low. Deep submm observations would 
be of value to better determine the dust content of SBS 0335-052. 




Fig. 17. — Dust-to-gas mass ratio as a function of Hubble type, assuming Xco = 4 x 10 20 cm~ 2 /(Kkm s _1 ). The 
horizontal line corresponds to M^uBt/Mu = 0.0073, and the broken lines to factor-of-two variations around this. Note 
that the upper bounds are all cases where the dust mass has been estimated, but only lower bounds are available for 
the gas mass because CO has not yet been observed. 



7.3. Dust-to-Gas Mass Ratio and Galaxy Type 

The dependence of the dust-to-gas mass ratio on galaxy type is explored in Figure [TTJ Spiral 
galaxies (Sab, Sb, Sbc, Sc, Scd, Sd) routinely have M^^/Mu « 0.007 to within a factor of two. 

We have dust mass estimates for 3 elliptical galaxies (NGC 0855, NGC 3265, and NGC 4125), 
but lack information on the gas content of NGC 0855 and NGC 4125. Detection of 21cm emission 
gives a lower bound on the gas mass in NGC 3265. If the molecular and ionized gas mass in 
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NGC 3265 is small compared to the HI mass, then this elliptical galaxy has a dust/gas ratio 
consistent with ~solar metallicity. 

The three SO galaxies in our samples have dust /gas ratios that are consistent with the range 
seen for normal spirals: NGC 3773 has M dust /M H < 0.01 and NGC 1482 and NGC 5866 have 
M dust /M H « 0.005. 
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Fig. 18. — The dotted histogram (shaded) shows the 17 galaxies with SCUBA data for which we can estimate the 
PAH index qpah- The solid histogram includes 44 additional galaxies without SCUBA fluxes for which (/pah can be 
estimated. Four galaxies for which qpah cannot be determined have been omitted (see text). 



7.4. PAH Abundance 

The set of 65 SINGS galaxies includes galaxies with strong PAH emission detected in IRAC 
bands 3 and 4, as well as cases where there is little evidence of nonstellar emission in these bands. 
For some galaxies the dust emission (relative to starlight) is so weak that qpah (the fraction of 
the total dust mass contributed by PAHs with Nq < 10 3 ) cannot be reliably determined from the 
broadband photometry. The DL07 models that are used here have 

WW. » 0.070 . (20) 
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With the existing uncertainties in IRAC calibration and foreground/background subtraction, we 
do not attempt to estimate (/pah unless 



(21) 



so that the PAH power in the IRAC 7.9/Lim band would be >50% of the stellar contribution if 
(/pah = 0.01. This implies no estimate of (/pah for four galaxies - NGC 1291, NGC 1316, NGC 4125, 
and NGC 4594 - all galaxies with relatively large bulge/disk ratios. Figure [18] shows the distribution 
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Fig. 19. — Histogram of global (/pah for SINGS galaxies with > 8.1, (a) with, and (b) without an AGN 
component. 



of (/pah for the remaining 61 galaxies for which we have both IRAC and MIPS data. The distribution 
is broad - 7/61 of the galaxies have (/pah < 0.75%, while 30/61 of the galaxies have (/pah > 3.25%. 



Roche et al.l (]199ll ) noted that the 8-13/mi spectra of AGN were often deficient in the PAH 



emission features commonly seen in H II region galaxies. The spectrum of PAH emission around 
AGNs has since been found to often exhibit unusu al band ratios, wi th some AGN systems show- 
ing extremely low values of L(7.7/mi)/L(11.3//m) (jSmith et al.l 120071 ). The suppression of 7.7/im 
emission in AGNs may be attributable to destruction of small PAHs by, e.g., X-rays from the AGN 
(jVoitl Il992l ). Figure [19] shows the distribution of global values of (/pah for SINGS galaxies with 
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low-lumin osity AGN (Seyfe r t or L INER), or (b) H II nuclei [the AGN/H II nuclear classifi- 



cation is from iMoustakas et al.l (|2007l )]. Only galaxies with characteristic metallicities Aq, > 8.1 
are included. The galaxies with AGN tend to have higher values of (/pah than the H II galaxies: 
the global SEDs show no evidence for PAH suppression. This is presumably because the SINGS 
sample excludes powerful AGN, and low-luminosity AGN have little effect on the global SED. Note 
that the 5 galaxies in Figure [T9b with (/pah < 2% include two irregulars (NGC 1705, Mrk 33), one 
elliptical (NGC 0855), and a dwarf starburst (Tol 89). 
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Fig. 20. — PAH index (percentage of dust mass contributed by PAHs with iVc < 10 3 C atoms) vs. galaxy metallicity 
(see text). Low metallicity galaxies always have low PAH index qpah ■ Filled circles are SINGS-SCUBA galaxies; 
diamonds are SINGS galaxies lacking submm data. 



Figure [20] shows (/pah versus metallicity for the 61 galaxies in the present study for which (/pah 
has been determined. We are able to estimate (/pah for galaxies spanning a range of at least 0.9 
dex in metallicity, from A Q = 7.68 ± 0.03 (Ho II) to 8.60 ± 0.06 (NGC 3351). [A Q = 7.54 ± 0.34 
of Ho I is nominally lower than for Ho II, but is uncertain.] A clear separation is seen: all nine of 
the galaxies with < 8.1 have (/pah < 1.9%, with median (/pah = 1.0%. When ^ 8.1 , higher 
PAH abundances are seen, with median (/pah = 3.55% for the 52 galaxies with Aq > 8.1, although 
note that only three of the 52 galaxies with Aq, > 8.1 have (/pah < 1-0% (the median for Ao < 8.1). 
PAH emission strength - as a fraction of total IR emission - therefore appears to always be low for 

< 8.1. 
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The weakness of PAH emission from low-metallicity galaxies was first noted by iRoche et al 



19911) using ground - based spec troscopy, and fu rther investigated with ISO data by iBoselli et al 



Garnett 



199il.lsturm et ali fcoooh . and lMaddenl (|200oh . The SMC, with A Q = 8.0 fcurt k Dufourlll998l: 



19991 ) falls just below the apparent threshold Aq = 8.1 seen in Figure [20l iLi &: Draine 



(I2002r) concluded that the SE P of the SMC Bar was consistent with the "SMC" dust model of 
Weingartner &: Draine ( 2001a ). with a PAH index q pah of only 0.1% E (see Tabled]). ISO observed 
PAH emission from molecular cloud SMC Bl No. 1 (jReach et al.ll2000l ) outside the SMC bar; Li & 
Draine showed that the observed emission could be reproduced with a model with only 3% of the 
SMC carbon in PAHs - corresponding to (/pah ~ 0.8%. Thus the SMC appears to have (/pah < 1%- 



Hunt et al.l (|2005l ) observed PAH emission to be weak in a number of low metallicity blue 
compact dwarf galaxies. In a photometric study of 34 galaxies spanning two decades in metallicity, 
Engelbracht et al.1 ()2005l ) concluded that there was a sharp difference in the ratio of 8/xm emission 
from PAHs to 24//m em i ssion from warm dust for metallicities above and below a threshold value 
Ao ~ 8. 20 iHogg et al. ( 2005 ) compared PAH emission with starlight for a sample of 313 SDSS 
galaxies, and noted that low-luminosity (presumably low metallicity) st ar-forming galaxies t ended 
to have low ratios of PAH emission to starlight. Using ISOCAM data, iMadden et al.l (|2006l ) con- 
cluded that PAH abundance, relative to larger grains, was positively correlated with metallicity, 
and suppressed in systems with hard radiatio n fields (i.e . , high [Ne IIl]/[Ne II] ratio). In a spectro- 
scopic study of blue compact dwarf galaxies, IWu et al.l (|2006l ) found that the equivalent width of 
the PAH 6.2 and 11.2/im features appeared to be suppressed for Z/Z & < 0.2, or Aq < 7.9. 



The threshold Ao = 8 .1 found here is c lose to the value 8.2 found by Engelbracht et al.1 (j2005l ). 
and the value 8.0 found bv lWu et al.l (J2006J). The reason for the small value of (?pah when Ao < 8.1 
is not clear at this time. A number of possibilities exist: 



Low Ao might imply more rapid destruction of PAHs by UV radiation in an interstellar 
medium with reduced shielding by dust. 

Low Ao might imply more effective destruction of PAHs by thermal sputtering in shock- 
heated gas that cools slowly because of the reduced metallicity. 



Galaxies with low O/H have lower values of C/O (jHenrv fc Worthevlll999l ); the reduced PAH 
abundance in these galaxies might be due to a deficiency of PAH-producing carbon stars and 
C-rich planetary nebulae. 



7 This result is controversial: for the SMC Bar, iBot et al.l ( 2004 ) have argued that the PAH abundance as a 
fraction of the dust mass is similar to the MW value - which we estimate to correspond to gpAH ~ 5%. Observations 
of the SMC with Spitzer should soon clarify this issue. 

8 Note that the oxygen abundances used by Engelbracht et al. (2005) were based on a heterogeneous compilation 
of measurements from the literature, whe reas the abundances in our study have been derived self-consistently, and 
placed on a common abundance scale (see iMoustakas et al.ll2007l . in prep., for details) 
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PAH formation and growth in the interstellar medium would be suppressed by low gas-phase 
C abundances. 

We also call atte ntion to the relatively low value of qpAn ~ 1-3% for Mrk 33, which with 



Aq = 8.30±0.10 (jMoustakas et al.ll2007l ) is well above the threshold A Q « 8.1. The relatively 
low value of qpAn in Mrk 33 may be a consequence of the vigorous star formation in Mrk 33 
which - via the combined effects of hard UV and supernova blast waves - may be enhancing 
the rate of destruction of PAHs, thus acting to lower the steady-state PAH abundance. 

• The strong 7.6/^m emission feature requires free- flying PAHs with Nc < 10 3 C atoms (DL07), 
and emission can therefore be suppressed if small PAHs coagulate with larger grains. Such 
coagulation can take place in molecular clouds. (How this might correlate with metallicity is 
uncertain, as the metallicity-dependence of molecular cloud properties is not understood.) 

While we can speculate about various effects that may be important, our limited understanding 
of the dynamics of the interstellar medium in other galaxies, together with our limited knowledge 
regarding many of the processes acting to form and destroy PAHs in the interstellar medium, 
precludes definite explanations of the observed galaxy-to-galaxy variations in (/pah- 



8. Starlight Properties for 65 SINGS Galaxies 
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Fig. 21. — C/ m in and (U) for 65 galaxies (see text). The shaded portion of the histogram consists of the 17 galaxies 
for which SCUBA global fluxes constrain the dust models. 
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The distribution of U m { n values is shown in Figure [2T1 The shaded portion of the histogram is 
contributed by the 17 galaxies for which SCUBA data were available and used in the model-fitting 
with fixed C/ max = 10 6 ; for these galaxies, there was no lower limit imposed on U m \ a . Two of the 
17 galaxies have U m i n = 2.0, and the other 15 have £7 m i n > 2.5. The remaining 48 galaxies lack 
submillimeter data, and for these the restricted fitting procedure allows only U m \ n > 0.7. This 
lower limit on U m \ n results in a small "pile-up" of five galaxies at U Ui \ a = 0.7 in FigureET]- possibly 
some of these should have smaller values of U m i n . What is needed, of course, are submillimeter 
observations so that the restricted fitting procedure becomes unnecessary. 

The form of the starlight distribution function (eq. [6]) is based on the expectation that the bulk 
of the IR luminosity would be radiated by grains in a diffuse interstellar medium with a more-or-less 
uniform starlight background, with U ~ U m ™ ■ This assumption appears to be validated by the 
success of the model in reproducing the observed dust emission extending from 5.7/xm to 850/um. 
It is interesting that the median U m i n = 1.5 is only slightly more intense that our local interstellar 
radiation field (U = 1), suggesting that interstellar processes (star formation, absorption of starlight 
by dust) regulate the diffuse starlight intensity to have U ~ 1.5 to within a factor of two. 

The distribution of mean starlight intensities is shown in Figure 121b . As before, the galaxies 
with SCUBA data are shaded. The galaxies with SCUBA data and those without SCUBA data 
both show a tendency to have (U) ~ 2.5, although about 12% of the sample has (U) > 10. 

The contribution of high intensity regions to the dust heating is shown in Figure [22] In many 
of the sample galaxies, the observed 24/im emission appears to require that a significant fraction 
of the dust heating take place in regions with high starlight intensities. Figure [22] shows that the 



7o 



median galaxy has ~8% of the total IR power originating from regions with U > 10 , and 
of the IR power originating from regions with U > 10 3 . In most cases these high-intensity regions 
are probably PDRs bounding dense H II regions. For example, the dust just outside t he ionization 



front in the Orion Bar PDR is heated by starlight with U ~ 3 x 10 4 ( Allers et al.ll2005l ). and similar 
regions should be present in other star- forming galaxies. For steady star formation with the IMF 
observed in the Milky Way, perhaps ~10% of the bolometric luminosity is contributed by the 
massive (M > 2OM ) stars that tend to spend most of their luminous existences close to their 
natal clouds, photodissociating and photoionizing the gas, and heating the dust. Therefore the 
median galaxies in Figure [22] do not seem surprising. 

The galaxies with the highest values of f(U > 10 3 ) are Mrk 33 (Im), Tolo 89 (Im), NGC 5408 
(Im), NGC 2798 (Sa), NGC 3049 (Sab), and NGC7552 (Sab). All of these galaxies appear to be 
forming stars at a high rate per unit area (see, e.g., discussion of Mrk 33 in §5.6.21) . and all have 
high ratios of ( vF u ) 24^ m /{yFv)i60nmi ranging from 0.52 for Tolo 89 to 1.47 for Mrk 33, indicative 
of a substantial fraction of the dust power originating in regions where the grains are hot enough 
to radiate at 24/im. 



DL07 showed that, for their dust models, the fraction of the starlight heating contributed by 
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(b) f(L d ;U > lxlO 2 ) 
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(c) f(L d ;U > lxlO 3 ) 
median = 0.062 
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Fig. 22. — Fraction of the dust luminosity originating in regions with (a) U 
(see text). Galaxies with / < 0.01 have been grouped in bin at 0.01. 



(b) Z7 > 10 2 , and (c) U > 10 3 



high intensity regions could be related to the observed quantity 



■21 



■0.14P 79 .^- aM( ^ nS) 
(vF v ) 71 + (vF v ) im 



(22) 



where F™ s is the nonstellar contribution to the flux density. The term 0.14(i/F" s )7.9 in the numerator 
of eq. (I22p is intended to approximate the contribution to the 24/im flux arising from single-photon 
heating of PAHs, so that the numerator of eq. (|22p is the 24^m flux contributed by larger grains, 
which will only radiate at 24//m when exposed to high intensity starlight. Figure [23] shows the 
estimated fraction of the dust power originating in regions with U > 10 2 vs. the observed value of 
the quantity (P24 — O.I4P7.9). The galaxies fall in an approximately 1-dimensional sequence, and 
the fraction of the starlight absorption taking place in high-intensity PDRs can be estimated from 
the observed (F" s )7.g, (i ? ° s )24, {F u )fi, and (F u )iqq using the DL07 fitting function 

f{L d ; U > 10 2 ) « 1.05(P 24 - 0.14P 7 . 9 - 0.035) 75 . (23) 
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9. Discussion 

9.1. What About the Dust in Dark Clouds? 

The adopted "model" for the distribution of starlight intensities makes no allowance for the 
possibilility of dust in dark clouds, since it assumes that there is no dust heated by starlight inten- 
sities below the "diffuse ISM" value C/ m in- While such dust may not be important as a contributor 
to the total bolometric luminosity, it would not be surprising if it contained an appreciable fraction 
of the total dust mass, and it might have been expected to contribute an appreciable fraction of 
the emission at 850/iin. 

However, we have seen that dust with the physical properties adopted by DL07 appears to 
reproduce the observed SEDs of the 17 SINGS-SCUBA galaxies without a "cold dust" component 
- indeed, when we fix U max = 10 6 , the best-fit models for the 17 SINGS-SCUBA galaxies all have 
U ■ > 2 

^min _ *>■ 

As discussed by DL07, the dust emissivities employed here reproduce the FIR and submm 
emission from our local diffuse ISM at galactic latitudes |6| > 7° as measured by the FIRAS instru- 
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ment on COBE (jFinkbeiner et al.lll999l ). To the extent that our models reproduce the observed 
global FIR and submm emission from other galaxies, it suggests that the dust in dark clouds makes 
only a minor contribution to the emission. Suppose that a fraction /cold of the dust mass is actually 
in dark clouds, with dust temperature T co \d, and the remainder (1 — / co id) is in diffuse clouds, all 
at a single dust temperature T^q. Then the fraction p co \d of the total emission contributed by the 
dark cloud dust at wavelength A is just 

/ cold /(e Xc - 1) 

Pc ° ld(A) = /coid/(e- - 1) + (1 - /coid)/(e- - 1) ' (24) 
where x c = he/ \kT co \& and Xd = he/ XkT^g. If we take = 20 K and T co id = 12 K, for example, 
then at 850/im we have 

p /cold (25) 

/cold + 2.33(1 — /cold) 

Hence if, hypothetical^/, 50% of the dust mass were in the cold clouds (/ co id = 0.5), they would 
contribute only 30% of the total emission at 850/im; if this hypothetical cold component were left 
out, the model would fall short at 850/zm by only 30%, comparable to the existing photometric 
uncertainties. 

Because the model does not include a cold dust component, it would not be surprising if the 
resulting estimates of the dust mass Md ust were low. However, we have seen that if we adopt 
Aco = 4 x 10 20 cm _2 (Kkm s -1 )" 1 , we obtain dust-to-gas mass ratios that appear to account for 
all of the refractory elements in the ISM, as shown in Figure [16j If in fact we have substantially 
underestimated the dust mass in some of these galaxies, then our adopted gas masses would have 
to also be increased, by using an even larger value of Xqo- 

Note that while much of the gas in the Milky Way is molecular, much of this molecular gas is 
actually in regions with visual extinction Ay < 1, so that it is not "dark", and the dust grains are 
heated by diffuse starlight. In dark clouds that are star-forming, such as Orion OMC-1 or M17, 
much of the dust is actually warm, heated by the embedded stars. It appears most likely that 
the fraction of the dust mass that is genuinely "cold" is probably relatively small (less than, say, 
30%). Although the model omits a "cold dust" component, it appears to recover the majority of 
the global dust mass in the SINGS galaxies. 



9.2. Comparison with Graphical Methods 

The model-fitting carried out here has been done by an automated fitting procedure, that 
includes using interpolated models with values of (/pah intermediate between the seven (/pah values 
in the Milky Way dust models examined by DL07. DL07 presented a straightforward graphical 
procedure for estimating (/pah> dust masses and starlight intensities from IRAC and MIPS pho- 
tometry, and applied it to 5 galaxies as examples. The results of that graphical procedure are 
compared with the present results in Table [71 For two of the galaxies (Mrk 33 and NGC 3521), the 
SED model fitting made use of SCUBA data, which is not used in the DL07 graphical procedure. 
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Table [7] shows generally good agreement between the graphical method and the present method 
for estimation of (/pah ; -^dustj an d (U). The dust masses agree to within a factor of two in all 
cases. The values of (/pah are in generally good agreement except for Mrk 33, where (/pah = 3.2% 
estimated from the graphical method is considerably larger than (/pah = 1-3% from the SED 
fitting procedure used in this paper. The reason for this discrepancy is straightforward: the DL07 
graphical fitting procedure advises the user to estimate (/pah by seeking a model that reproduces 
the 7.9fim emission with the smallest possible value of 7, the fraction of the dust mass exposed to 
starlight with U > U m i n . While this approach is appropriate for most galaxies - which generally 
have 7 < 0.03 - Mrk 33 has an extreme spectrum, with a very large fraction of the dust exposed 
to intense radiation fields (see the location of Mrk33 in Fig. [22]) . The SED modeling in the present 
paper estimates 7 = 0.12 for Mrk 33. In fact, examination of the graphs shown in DL07 (Figs. 20 
and 21 of that paper) reveal that a solution with 7 ~ 0.1 does provide a better fit to the observed 
flux ratios. The DL07 graphical fitting procedure should therefore be employed with caution for 
galaxies with P24 = {vF 1 £ s )2±/\{vF v )i\ + (^-^v)i6o] ^ 0.2, indicative of substantial amounts of IR 
power radiated by dust in intense radiation fields. 



9.3. Dust Composition 



The graphite component of the dust model used here pro duces a broad emis sion feature near 
33/im when the dust is illuminated by starlight with U > 10 4 ( Draine &: Lil 120071 ) . This feature is 
a consequence of the optical properties assumed for graphite. 



Visual inspection of the 25-35/xm spectra presented by ISmith et al.l (|2007l ) suggest that the 
33/xm feature is not present, at least not with the strength that would be expected in some of 
our global models. Some of the central spectra (e.g., NGC 4321) have a broad (FWHM ~ 1.7/xm) 



Table 7. Comparison with Results of Graphical Method 
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1.1% 


0.7% 
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5.04 


2.1 


2.8 



Graphical method from DL07 

Results from fitting MW models with U max = 10 6 and U m \n > 1. 



55 



emission feature near 33/xm, but it is weak0 and narrower than the very broad (FWHM ~ 20//m) 
peak found in our models. A similar ~ lfim wide platea u near 33-34/im was als o seen in the 



ISO spectra of bright starburst galaxies ([Sturm et al.l |2000| ) and planetary nebulae ([Waters et al 



19981 ). and tentatively attributed to crystalline olivine. If future observations find that the graphite 
feature is not present with the strength predicted by our model calculations, what would be the 
implications? First of all, it must be remembered that the present calculations are based on two 
highly idealized assumptions: (1) the a > 100 A carbonaceous grains are spherical and composed of 
single crystals of graphite, and (2 ) the dielectric tensor in the IR conforms to a simple free-electron 
model iDraine &: Led (jl984[ . 119871 ) . The 33/im feature is a consequence of the dielectric function 
for E || c, which is directly related to the electrical conductivity parallel to the c-axis (i.e., normal 
to the graphite basal plane). This conductivity is modelled by a free electron gas with a plasma 
frequency lu p = 1.53 x 10 14 s _1 , and a damping time r = 1.4 x 10 _14 s. The simple free-electron 
model is not expected to accurately describe the detailed wavelength dependence of this feature 
(located at to = 27rc/33 / um = 5.7 x 10 13 s~ 1 ). Furthermore, even if the model were accurate for 
ideal graphite, it seems likely that impurities and defects - which might vary from grain to grain 
in a single region, or from region to region in the galaxy - would have the effect of causing oo p and 
r to vary over the population of emitting grains (and perhaps even within a single grain) resulting 
in further broadening of this feature. Therefore, failure to observe this feature would probably not 
rule out graphite as a grain constituent - it would only rule out monocrys talline graphit e grain s 
with a dielectric function resembling the simple free-electron model used bv lDraine &; Led ([1984]). 



Of course, it would not be surprising if the carbonaceous component of the interstellar grain 
population consisted of something quite different from monocrystalline graphite. "Amorphous 
carbon" , for instance, refers to a class of materials with widely var ying properties: for example , 
the amorphous carbo n materials BE, ACAR, and ACH2 studied by IColangeli et al.l (|l993l . Il995l ): 



Menne 



la et al 



( 19951) ar e all insulators, with bandgaps of 0.15 eV, 0.52 eV, and 1.32 eV respec 



tively (jZubko et al.lll996l ). These carbon solids, each having different properties in the FIR, would 
not show the 33^m feature. It is quite possible that the carbonaceous material in the larger 
(a > 0.02/um) interstellar grains may be better approximated by some form of amorphous carbon. 
However, we have seen that the SINGS photometry appears to also be compatible with graphite. 



9.4. The H 2 /CO Conversion Factor X C o 

To estimate the dust-to-gas mass ratio in §7.11 we need the total hydrogen mass Mh, including 
H2. Because H2 resides primarily in its lowest two rotational states (J = and 1), which do not 
radiate, the CO J = 1 — > line is routinely employed to trace H2. It is usually assumed that 
the H2 column density is proportional to the CO 1-0 surface brightness (see equation I14p . with an 



9 The PAHFIT spectrum decomposition program ( Smith et al. 
these galaxies. 



20071 ) assigns an equivalent width ~ 0.2/im for 
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empirical conversion factor Xqo that is based on some method of estimating the total gas mass. 
For the Milky Way, virial mas s estimates based on gia nt molecular cloud linewidths give Xqo = 
3.0 x 10 20 cm" ^(Kkms -1 )" 1 (jYoung fc Scovilld Il99lh . Studies using 7-ray e mission as a m ass 
trace r in the Milky Way fou nd Xqq / [ cm" 2 ( K k m s" 1 )^ 1 ] = (1.9 ± 0.2) x 10 20 ([strong & Mattox 



1996 ) and (1 . 56±0 .05) x 10 20 ( Hunter et al.lll9971 ). Using IR emission from dust as a mass indicator, 



Dame et alJ (|200lh found X CO = (1.8 ± 0.3) x 10 20 cm- 2 (Kkm s" 1 )" 1 for gas at 5° < \b\ < 30°. 



Because (1) CO 1-0 emission is often optically thick, (2) CO may be subthermally excited in the 
outer layers of molecular clouds, and (3) photodissociation of CO results in low CO/H2 ratios in 
the outer layers of molecular clouds, there is no reason to expect a single value of Xqq to apply to 
all galaxies, or to all regions within a galaxy. 



For low-luminosity spiral galaxies in the Virgo cluster, iBoselli et alJ (12002 ) estimated values 



of Xqq ranging from ~ 10 cm (Kkm s 



for giant spirals to ~ 10 cm (Kkm s 



for 

dwarf irregulars. B ased on an e xtens ive study of giant molecular clouds in a number of galaxies 
in the Local Group. iBlitz et al.l (120061 ) conclude that Xqo = 4 x 10 20 cm -2 (Kkm s -1 ) -1 is a good 
value to use for most galaxies if the individual GMCs are assumed to be virialized. 

The present study has demonstrated that we can use a physical dust model to reproduce 
observed IR and FIR emission from nearby galaxies, allowing us to estimate the dust mass in these 
galaxies. Because we have a-priori expectations for the value of the dust-to-gas mass ratio, our 
determination of M^nst therefore allows us to estimate the total gas mass, and thereby estimate 
Xqo for these galaxies. 

Figure [TBI explores the relationship between dust-to-gas mass ratio and metallicity. Observed 
values of Mj us t/MH are plotted versus metallicity in Figures [T5r and b, showing only galaxies where 
both H I and CO have been detected, or where the upper limit on either H I or CO shows it to 
be of minor importance compared to the other gas component. Equation (|13f) is plotted in Figure 
[T5a.b. where we assume [Ao]mw = 8.59. The figure includes 12 galaxies for which SCUBA data 
constrains the dust models - these should be the most reliable dust masses. The figure includes 
an additional 20 galaxies lacking submm data. In Figure [T5h the H2 mass is estimated using 
Xqq = 2 x 10 20 cm~ 2 (Kkm s" 1 ) -1 , whereas Figure [T5b uses Xqq = 4 x 10 20 cm~ 2 (Kkm s" 



From comparison of Figures [T5a and b, we see that Xqo = 4 x 10 20 cm~ 2 (Kkm s" 1 ) -1 results 
in dust-to-gas mass ratios for the SINGS-SCUBA galaxies that are consistent with the amount of 
material actually available to form silicate and carbonaceous grains, when one considers that the 
dust mass estimates are uncertain at the factor-of-two level, and the metallicity estimates are also 
uncertain. There are some cases with mass determinations (not upper limits) that fall well above 
eq. (fl~3l) . which was expected to define an upper envelope for M^t/Mu. For all cases where submm 
data were used in the fitting, the dust mass estimates were expected to have better than factor of 
two accuracy; if we adopt Xqo = 4 x 10 20 cm _2 (Kkm s^ 1 ) -1 , none of the galaxies with SCUBA 
data exceed eq. (|13p by more than a factor 2. Conversely, Xqo = 2 x 10 20 cm _2 ( K km s" 1 )" 1 would 
imply dust masses for NGC 1482 and NGC 2976 that would be uncomfortably large given their 
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estimated oxygen abundances. The one galaxy in Figure 115b that falls well above the expected 
envelope is NGC 4594, the Sombrero galaxy. The interstellar medium in this unusual system merits 
further study. 

The histograms in Figures [T5b.d show the distribution of the dust-to-gas mass ratio, normalized 
to the characteristic oxygen abundance of the Milky Way. Ideally, the distribution would be narrow 
and peaking at ~0.01, possibly with a tail to low values of (Md us t/MH)/[(0/H)/(0/H) MW ] repre- 
senting galaxies where the grains do not incorporate a large fraction of the refractory elements. For 
Xqo = 2x 10 20 cm- 2 (Kkm s -1 )" 1 , the median is at 0.016. For X co = 4 x 10 20 cm" 2 ( K km s" 1 )" 1 , 
the median shifts to 0.012 - approximately the value expected if grains incorporate most of the 
refractory elements (as is the case in the Milky Way). 

The a bove consideration s suggest that the value Xqq ~ 4 x 10 20 cm~ 2 (Kkm s" 1 )" 1 recom- 



mended by iBlitz et al,l (|2006l ) for Local Group galaxies is also appropriate for the SINGS galaxy 
sample. This valu e of Xco is a f actor of ~2 larger than the value (1.8±0.3) x 10 20 cm _2 ( K km s -1 ) -1 
found recently by I Dame et al.l (|200ll ) for the local Milky Way. 



10. Summary 

The principal results of this paper are as follows: 

1. A physical dust mo del consisting of PAHs, carbonaceous grains, and amorphous silicate grains 
(jPraine Lill2007h heated by starlight, with a distribution of starlight intensities, successfully 



reproduces the IR and submm emission for a sample of 17 galaxies with IRAC and MIPS 
photometry plus SCUBA photometry. In particular, the SCUBA photometry does not require 
an additional population of cold dust grains, or grains with enhanced submm emissivities. 
Cold dust grains in dark clouds appear to contain only a minor fraction (< 50%) of the total 
dust mass. 

2. The dust model allows the mean starlight intensity scale factor (U) and the total dust mass 
^dust to be estimated for each galaxy. For the 17 galaxies with SCUBA data, (U) ranges 
from 2.2 to 16, with a median of 4.3 (see Fig. [7]). For the 12 galaxies with SCUBA data 
and CO and H I observations, and assuming Xqo = 4 x 10 20 cm~ 2 ( K km s" 1 )" 1 , the median 
M dust /M u = 0.0052 (see Fig. gj) - this is similar to M dust /M H ~ 0.0073 estimated for the 
Milky Way based on observed depletions (see Table [2]), or the value ~0.010 estimated from 
dust models that reproduce the observed Milky Way extinction (see Table [3]) . 

3. When submm observations are unavailable, the mass of cool dust in a galaxy is poorly con- 
strained. However, we use the model fits for the 17 galaxies where SCUBA data are available 
to test a "restricted" fitting procedure that can be applied when submm observations are 
unavailable. For the 17 galaxies where global SCUBA data are available, the restricted fitting 
procedure recovered the dust mass to within a factor 1.5 in 11/17 cases (see Fig. \12\i . and to 
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within a factor 2.2 in every case. This restricted fitting procedure can therefore be used to 
estimate dust masses and starlight intensities in galaxies with dust emission detected in the 
three MIPS bands. 

4. For the 20 galaxies which lack SCUBA fluxes but for which the gas mass is known, the median 
M dust /M H « 0.0088 (see Fig. Q3J. 

5. By comparing dust-to-gas mass ratios computed with different assumed values of the CO to 
molecular mass conversion factor Xqo, we find that the dust-to-gas mass ratio normalized to 
galaxy metallicity gives dust masses consistent with available metal abundances for Xqo ~ 
4 x 10 20 cm- 2 (Kkm s -1 )" 1 (see Fig. EE]). 

6. All 49 galaxies in Figure [16] with Ao > 8.2 appear to have dust masses consistent (to within 
a factor ~3) with eq. (fT3l) . corresponding to the total mass of dust in the WD01 models for 
Milky Way dust (see Table [3]) scaled by metallicity. This is close to the mass that will be 
formed if silicates are formed from the bulk of the interstellar Mg, Si, and Fe, and ~50% of 
the interstellar carbon goes into carbonaceous grains (see Table [2]). 

All 9 of the galaxies in Fig. [16] with metallicities Ao < 8.1 appear to have global dust-to-gas 
ratios that fall below eq. (|13p . in some cases by large factors: NGC 2915 appears to be at least 
a factor 40 below eq. (fl~3j) . However, the dust-to-gas ratio in the regions where dust emission 
is detected appears to be close to eq. (fl~3j) . The low apparent global dust-to-gas ratios are 
found in dwarf galaxies with extended H I envelopes. When allowance is made for the lower 
metallicity and lower starlight intensity in the envelopes, it appears that the dust-to-gas mass 
ratios may be consistent with eq. (13). 

7. We characterize the PAH abundance by a quantitative "PAH index" (/pah, defined to be the 
fraction of the total dust mass contributed by PAHs containing less than 10 3 C atoms. The 
61 galaxies in this sample for which we are able to estimate (/pah span the range of PAH 
index from 0.4 to 4.7%, with a median (/pah = 3.4% (see Fig. [T8l) . 

8. The PAH index (/pah is strongly correlated with metallicity (see Fig. [20]) . The nine galaxies 
with Ao < 8.1 have (/pah < 2%, and median (/pah = 1-0%, whereas the 52 galaxies with 
Ao > 8.1 have median (/pah — 3.55%. Metallicity Ao = 8.1 appears to mark a transition in 
the composition of interstellar dust in galaxies, from low- PAH to high- PAH. Possible reasons 
for the observed variation in (/pah are discussed in £ 17.41 

9. Except for three starbursting galaxies (Mrk 33, Tolo 89, NGC 3049), the infrared emission 
from dust is dominated by dust grains exposed to starlight intensities U ~ U m \ a , interpreted 
here as dust in diffuse regions heated by the general diffuse starlight background (see Fig. 
[22h). Half of the galaxies have more than 88% of the total dust power provided by this 
component. 

10. The SEDs for the star-forming galaxies in the sample require that a significant fraction of 
the dust luminosity be produced by warm dust in regions of intense starlight. Half of the 



59 



galaxies in the sample have more than 8% of the total dust luminosity originating in regions 
with U > 10 2 (see Fig. [22]) . This appears to be consistent with the expected heating of dust 
in PDRs near massive stars. 

We thank S. Madden and the anonymous referee for helpful comments. Support for this work, 
part of the Spitzer Space Telescope Legacy Science Program, was provided by NASA through the 
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Spectral Energy Distributions for 48 Galaxies Lacking Submm Data 
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